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We use QCD to compute the cross section for high-energy coherent production of a dijet (treated 
as a qq moving at high relative transverse momentum, Kt) from a nucleon and a nuclear target. 
The direct evaluation of the relevant Feynman diagrams shows that, in the target rest frame, the 
space-time evolution of this reaction is dominated by the process in which the high Kt qq component 
(point like configuration) of the pion wave function is formed before reaching the target. This point 
like configuration then interacts through two gluon exchange with the target. In the approximation 
of keeping the leading order in powers of a s and in the leading log approximation in a 3 ln(fvt /^qcd)i 
the amplitudes for other processes are shown to be smaller by at least a power of a„ and/or powers 
of Sudakov-type form factors and the small probability, W2, to find a qq pair with no gluons at 
an average separation between constituents. Thus the high Kt component of pion wave function, 
including the contribution of Gegenbauer polynomials of rank n > 0, can be measured in principle at 
sufficiently large values of k\ . At large values of k\ , the resulting dominant amplitude is proportional 

to z(l— z)a s (kt )« t (In -p-) - ^ (z is the fraction light-cone (+) momentum carried by the quark in the 
final state, (3 is the coefficient in the running coupling constant) times the skewed gluon distribution 
of the target. For pion scattering by a nuclear target, this means that at fixed xn = 2nf/s (but 
k^ — > oo ) the nuclear process in which there is only a single interaction is the most important one to 
contribute to the reaction. Thus in this limit color transparency phenomena should occur-initial and 
final state interaction effects are absent for sufficiently large values of Kt- These findings are in accord 
with the recent experiment performed at FNAL. We also re-examine a potentially important nuclear 
multiple scattering correction which is positive and varies as the length of the nucleus divided by an 
extra factor of l/nf ■ The meaning of the signal obtained from the experimental measurement of pion 
diffraction into two jets is also critically examined and significant corrections are identified. We show 
also that for values of Kt achieved at fixed target energies, dijet production by the electromagnetic 
field of the nucleus leads to an insignificant correction which gets more important as Kt increases. 
We explain also that the same regularities are valid for photo-production of forward light quark 
dijets. 



I. INTRODUCTION 



The theory of strong interactions, QCD, contains many specific predictions regarding the space-time evolution of 
high energy coherent processes. A review describing the many interesting qualitative results that have been obtained 
in this rapidly developing field is provided in Ref. jl[ . The aim of this paper is to use a specific example of a completely 
calculable process to demonstrate the general properties of space-time evolution of hard exclusive processes in QCD. 
In particular, we consider a process in which a high momentum (~ 500 GeV) pion undergoes a coherent interaction 
with a nucleus in such a way that the final state consists of two jets (JJ) (formed by a qq pair) moving at high 
transverse relative momentum greater than about 2 or 3 GeV. The process of two-jetproduction was first discussed 
for both photon and pion projectiles interacting with a nucleon target 0], and Ref. || introduced the possibility of 
using this process to probe the nuclear filtering of small color dipoles. Estimates of Ref. || found that for heavy 
nuclei, nuclear filtering causes the exclusive dijet production to decrease exponentially as k\ increases. Hence an 
overall increase of the total diffractive cross sections was suggested as a good signature of the nuclear filtering of small 
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size configurations 0. In || we presented the first application of QCD to the process of dijet production at large Kt 
by generalizing QCD factorization theorems, predicted a nuclear dependence which is qualitatively different from that 
suggested in || , and a Kt dependence which differs by a power of k 2 from Q and qualitatively from that discussed in 
H . We also argued that this process can be used to directly measure the behavior of the qq component of the pion's 
light cone wave function for large values of K t . 

If one wishes to describe hard diffractive processes, it is important to realize that the effective number of bare 
particles in the light cone wave function of the projectile depends strongly on the longitudinal distances involved. 
If these distances are small, and the process involves high momentum transfer, (as occurs in computing hadron 
electromagnetic form factors) the main contribution originates from the Fock component of the hadron wave function 
containing the minimal number of constituents. On the contrary, if the longitudinal distances are sufficiently large, 
the minimal Fock component (a qq pair in our case) will develop additional components such as qq pairs and gluons. 
Thus for those processes in which infinite longitudinal distances are involved, the number of partons in the light cone 
wave function of the projectile would always be infinite. For processes initiated by a spatially small colorless dipole, 
using the QCD factorization theorem allows one to trace the origin of wee partons (carrying a small fraction of the 
momenta) as arising from the space-time evolution of the projectile's minimal Fock component, and to include these 
effects in the (skewed) parton distribution of the target. 

We now discuss the basic physics of the dijet production process. The selection of the final state to be a qq pair 
plus the nuclear ground state causes the qq component of the pion wave function to dominate the reaction process. 
At very high beam momenta, the pion breaks up into a qq pair with large K t well before hitting the nucleus. The 
dominance of this starting point is verified in the present work. Note also the crucial feature that for values of x which 
are not very small, so that the leading twist approximation is valid for the small dipole - target interaction, a spatially 
small wave packet of quarks and gluons remains small as it moves through the target. This leads to a dominance of 
the effects of large transverse momenta, and allows the factorization of the hard physics from the soft physics. On the 
other hand, for very small values of x, the packet lives so long that it would expand to a normal hadronic size causing 
the initial state interaction to become similar to the soft one. Moreover, a rapid increase of PQCD amplitudes with 
energy leads to the violation of QCD evolution equation at rather small values of x, and to the disappearance of the 
characteristic physics of the interaction of a small dipole. Thus QCD predicts different calculable dependencies of the 
cross section of the diffractive dijet production on atomic number, on n t and z in different regions of the \iik 2 /Aq CD , 
In 1/x plane. 

Thus two large parameters s and Kt 2 are present, and this feature will enable us to demonstrate the dominance of 
Feynman diagrams of a very few specific topologies, and to evaluate them. The result of calculations can be represented 
in the form of a generalized QCD factorization theorem, valid for the set of Feynman diagrams corresponding to the 
leading power of s and the minimal power of 1/k 2 and a s at fixed values of ot s In k\ Aq CD . Our calculation relies 
heavily on the well-known theoretical observation that amplitudes of many high energy processes (such as that in 
the parton model, multi-peripheral processes and those involving the Pomeron) are dominated by ladder diagrams 
HQ. This property has been proved using the approximation of including the terms of lowest order a s and all 
powers in a s InQ 2 /Aq CD |||| and/or a s lnl/a; |p|,^0[. The dominance of ladder diagrams makes it possible to 
absorb the effects of the leading terms in Iylk 2 /Aq CD and Inl/x into the dipole-target interaction, and/or into the 
target's skewed parton distributions. Furthermore, one may classify and analyze those diagrams of leading power of 
a s which are relevant for the pion transition into two jets. This leads to a selection rule: the i-channel exchanges 
with vacuum quantum numbers (positive charge parity) should dominate. Thus we will calculate an amplitude which 
is symmetric under the transposition s <-> u. Negative charge parity contributions (such as effects of the odderon) are 
neglected. There is another group of correction terms of the form a s In k 2 / Aq C d which arises from radiative effects in 
the pion wave function. The ladder structure of the dominant diagrams makes it possible to include these terms into 
QCD evolution of pion wave function and structure functions of the target will not change the structure of formulas 
deduced by considering skeleton diagrams calculated within the approximation of keeping leading powers of a s ; their 
only influence is to introduce the effects of evolution in k\ into the relevant parton distributions and into high k t 
behavior of the pion wave function. Indeed, our theoretical analysis heavily relies upon specific properties of skewed 
parton distributions and of minimal Fock component of pion wave function. Note also that, for small values of x, 
and large values of Q 2 , the skewed parton distribution of a target nucleon or nucleus is calculable in QCD using the 
appropriate evolution equation and initial diagonal parton densities fllj]. To simplify the calculations and especially 
the separation of the scales, we use an axial light-cone gauge, which reduces to the A- = gauge, in the target 
rest frame. This gives a high n t behaviour of fermion propagator and hard gluon exchange amplitude which have no 
infrared singularities. In this gauge, unphysical degrees of freedom are removed from the light cone pion wave function 
at least within the leading log approximation. Consequently, the separation of momentum scales can be easily made. 
On the contrary, in the gauge used in (4(| the separation of scales and therefore derivation of the QCD evolution 
equation, although correct, is complicated by the need to account for the cancellation of the infrared divergences. 



2 



We want to stress that in the calculation of amplitudes of hard diffractive processes in the gauge where fermion 
propagators are infrared divergent one, should first remove infrared divergences and only then consider partons in the 
non-perturbative wave function to be on-mass-shell. 

These technical considerations lead to some simple results for situations, such as ours, in which the momentum 
transfer to the target nucleus is very small (almost zero for forward scattering). In this case, the dominant source 
of high momentum must be the gluonic interactions between the pion's quark and anti-quark. This is also justified 
in the present work. Because is large, the quark and anti-quark must be at small separations-the virtual state of 
the pion is a point- like-configuration [}12| . But the coherent interactions of a color neutral point-like configuration are 
suppressed (at fixed xjsr, k| — > oo), for the processes which involve small transfers of momentum to the target, by the 
cancellation of gluonic emission from the quark and anti-quark p3|,p|,^2[ and /or from the qqg state. -see the discussion 
in the section II. D. 

Furthermore, the strength of the interaction with the target is proportional to the square of the transverse separation 
distance between the quark and anti-quark. Thus the interaction with the nucleus is very rare, and the pion is most 
likely to interact with only one nucleon. The result is that in this process the initial tt and the final qq pair do not 
get absorbed by the target, as would typically occur in a low momentum transfer process. Thus initial and final 
state interactions are suppressed and color transparency unambiguously follows. As the values of x are decreased, the 
qualitative physics changes gradually. The increase of the effective size of the small color-dipole leads to the increase 
of the influence of initial state interactions, and to a contribution of nuclear shadowing which enters at leading twist. 
For even smaller values of x, the leading twist approximation breaks down. 

Our treatment of the reaction process in terms of a separate wave function and interaction pieces provides a new 
example of how the QCD factorization theorem works for high energy processes involving two large variables: s, K t 
[ p2| . For this coherent process, the forward scattering amplitude is almost proportional to the number of nucleons, 
A, and the cross section varies as A 2 . The forward angular distribution is difficult to observe, so one integrates the 
angular distribution, and the A 2 variation becomes « Rii ^ +Bn oc 1+0 45^-2/3 ~ A 137 . Here Bm ~ A.hGeV~ 2 is the 
slope of the t dependence of the cross section of a hard diffractive process as determined by data for electroproduction 
of vector mesons. This very rapid variation represents a prediction of a very strong enhancement which occurs via the 
suppression of those interaction processes which usually reduce the cross section. Our interest in this curious process 
has been renewed recently by exciting experimental progress |l4| . 

Three key predictions of our paper || are confirmed by the E791 data: 

• The result from the E-791 experiment comparing Pt and C targets that the coherent cross section for small 
momentum transfer to the nucleus varies as ~ ^l.ssio.os^ j s c i ose to our predictions, see Section V. This 
variation is much stronger than seen in soft diffraction of pions by nuclei ~ ^4 ' 8 fl^ , |i~5|| , which is qualitatively 
different from the behavior ~ A 1 / 3 suggested in Ref. This A-dependence is somewhat more rapid than that 
predicted by color transparency theory for A — > 00. For moderately large values of A, small effects discussed 
in Sect. V tend to increase the A dependence. This may be understood as a result of the experimental trigger 
not excluding a small but calculable admixture of the effects of nuclear disintegration processes which lead to 
similar dependence of cross section on t. Furthermore, it is an unusual feature J5| of the present process that 
final state interactions of the point-like configuration tend to increase the A-dependence. Sect. V also contains 
a discussion of the changes of the A-dependence due to the effects of nuclear shadowing on the gluon density. 

• The dependence of the cross section oc z 2 (l — z) 2 on the fraction of momentum z carried by one of the jets is 
consistent with our prediction. 

• The cross section da/dn 2 falls as n~[ n with n — 10.2±0.4 (stat)±0.3 (sys),for n t > 1.25GeV and as n = 7.5±2.0 
for Kt > 1.8 GeV. This should be compared to the prediction of n — 8 p). For smaller values of n t soft QCD 



phenomena, such as production of qqg jets, should be important. See the discussion in Sect. Ill A 



The purpose of the present work is to rederive and confirm our earlier theoretical results with a more extensive 
analysis. The derivation of our leading term || directly from QCD by generalizing a QCD factorization theorem of 
Ref. jl| was presented in Ref. fl7[] , and this is explained more fully now. But here we go further by verifying the 
assumption that the point-like-configuration is indeed formed well before the projectile reaches the nucleus. In the 
derivation we shall explain that several different amplitudes, which seem to be of the same or lower order in a s as 
the leading one described above really are very small after proper account of the suppression of radiation collinear to 
pion momentum direction. 

We also update our study of the leading multiple-scattering correction, which is positive because the strength of 
the final state interaction decreases with decreasing size of the dipolc || , and we study the most important competing 
electromagnetic process. Some specific features of the experimental extraction of the coherent part of the cross section 
are also explained. Still another feature involves the soft interaction between the dipole and the target. This was 
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at first derived to be proportional to the gluon density of the nucleus |18| . However, there is a non-zero momentum 
transfer to the nucleus, so it is actually the skewed gluon density that enters. The skewedness of gluon distribution in 
the nuclear target leads to a small, calculable correction to the predicted A dependence |llj and absolute value 
which changes the detailed nature of our results but not the qualitative features. 

Our main results are summarized by the following formula, valid in the leading log approximation, for the differential 
cross section of diffractive dijet production by nuclei: 

da(Tr + A^2jet + A)(g t =Q) (1 + r) 2 ) ( f j2 d0 

du^ t = M2^(j d r *j e ^ K * ■ r *> 

/Ja{xi xi (3s k 2 ) \ 2 
d 2 k u [2xn(z, n) - exp(-i(r t • kit)x,r(2 - P, n) - exp(i(r t • kit)x,r(2 + (3, n)] \\k 2 ' ~ J ' ^ 

where by definition rt is the transverse distance between the pion's quark and anti-quark, x J^^y = 

J Gpr 4 expz(k t • r t )xi,Tr( z , r t) and x\Ga{xi, %2, Q 2 ) = d 2 l J^ o ^-d 2 l fA( - Xl,x ^i^' 1 ^ fc " (/a can be denoted as the un- 

integrated skewed nuclear gluon density), ki are four momenta of two exchanged gluons, (see Fig. 1), and k\ t — Kt — it. 
Here df(n 2 ) is the renormalization factor for the quark Green function St(k) in the hard regime where 

s f (k)=m, ( 2) 

for k 2 — k 2 . The quantity /?o is a complicated function involving the transverse momenta of the quarks within the pion, 
and in the region giving the dominant contribution (3q oc k 2 /s. The quantity r\ is the ratio of the real to imaginary 
part of the gg-target scattering amplitude. In Eq. ([!]), %(z,rt) includes both the non-perturbative qq component 
and its high momentum tail. This function therefore involves distances significantly smaller than average hadronic 
inter-quark distances. The actual distances involved in the largest contributions is one of subjects investigated here. 

Equation (|l|) is derived in two steps. First we demonstrate the dominance of the Feynman diagrams of Fig. 1 and 
then evaluate these diagrams. The factorization of the hard perturbative QCD part, related to the pion wave function, 
and qq pair arising from softer QCD, described by skewed parton distributions (the dominance of the diagrams of 
Fig. 1) is another form of the QCD factorization theorem derived in |l6| for diffractive vector meson production in 

A 2 

deep inelastic scattering DIS. The end point contribution - the Feynman mechanism -(z oc Q % D ) is suppressed as 

K t 

compared to the leading term by a set of factors: one power of -i^, the square of the Sudakov-type form factor, by a 

form factor W2 (which accounts for the very small probability to find a pion with q and q at average distances without 

a gluon field) and by the overlap integral with final state. A detailed analysis of the end point contribution will be 

subject of a separate publication. It follows from QCD factorization theorems that the amplitude of hard processes can 

be represented as the convolution of the non-perturbative pion wave function and hard amplitude T. The virtualities 

of all particles in the s, u cuts of the amplitude T are large. Virtualities of those seemingly on-mass-shell particles are 

3> Aq CD but <C k 2 . This is the condition which dictates the dominance of perturbative tail in the pion wave function 

at k 2 — > oo but fixed x = n 2 /v. The amplitudes having different topology are radiative corrections involving extra 

powers of a s . To elucidate the underlining physics we shall prove the dominance of the amplitude T± (see Fig. 1) 

by analyzing different contributions of many diagrams. A more general proof will be given elsewhere. In the leading 

o 2 i .... . 

In -nr — and In - approximations the dipole description can be used to simplify the above equation to the form: 



da(ir + A^2jet + A)(q t = 0) (l + ?7 2 ) 



dtdzd 2 n t 167r(27r) 



A - 2 —z—x 1 Ga{xi,x 2 ,Q ) 
• d f {K t ) 



2 



(3) 



where x-niz, Kt) = 47tCf as( -" t - ) -\/3f w z(l — z), A is the Laplacian in Kt space, z is the fraction light-cone (+) momentum 

carried by the quark in the final state, x\Ga{x\, x 2 , k 2 ) is the skewed gluon density of the nucleus, xi,x 2 are the 
fractions of target momentum carried by exchanged gluons 1 and 2, x\ ~ X2 — Ml^/s, X2 < x%, and r\ — ReF/ ImF 

(with F as the dipole- nucleon scattering amplitude). Note that the resulting k^ 8 dependence is a consequence of a 
kind of dimensional counting, as explained in Sect. II. 

It is necessary to discuss the kinematic and dynamic limitations of our analysis. We require high beam energies so 
that the point-like configuration remains small as it passes through the nucleus, and we also require that Kt be large 
enough so that the qq pair actually be in a point-like configuration. This situation corresponds to k 2 /s being held 

fixed for large values of k\. For the experiment of Ref. |l4|] K t ~ 2 GeV, and s = 1000 GeV 2 , so xn = ~ -008. 
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There is another kinematic limit in which k 2 is fixed and s goes to oo. At sufficiently small values of xpj, less than 
about 2m ^ Ra , the situation is very different because the qq dipole system is scattered by the collective gluon field of 
the nucleus. Nuclear modifications (enhancement) of the nuclear gluon density actually occurs at larger values of xn 
corresponding to xn ~ l/(2mjv t^n) ~ 0.1 (where r^N ~ 2 fm is the mean inter-nucleon distance in nuclei [^0] p2|). 
But for values such that xm < 2m ^ R the nuclear gluon field is expected to be shadowed, leading to a gradual 

disappearance of color transparency (at a fixed scale (fcf)). This is the onset of perturbative color opacity J5|, p3| , p4[ . 
At even smaller values of xjy a new phenomenon has been predicted - the violation of the QCD factorization theorem 
p5[ . Our present analysis is not concerned with this region of extremely small xn ■ Another interesting phenomenon 
is the possibility of probing the decomposition of quark distribution amplitude in terms of Gegenbauer polynomials 
at sufficiently large values of k 2 . 

Some general features of our analysis appear in several different Sections, so it is worthwhile to discuss these here. 
The calculations of several amplitudes are simplified by the use of a general theorem. In the leading order in 1/k 2 , the 
interaction of the qq occurs via the exchange of a two-gluon ladder with the target. It is important to note that the 
interaction of the qq pair with the target via the exchange of a larger number of gluons is suppressed by powers of 1/k 2 . 
The proof of this statement follows and heavily uses Ward identities. (To visualize the similarity with the situation 
considered in Ref. Jl6|], it is instructive to neglect the effects of the odderon contribution, which is in any case small. 
Accounting for symmetry of amplitude on the transposition: s <-> u gives the possibility to consider the amplitude for 
the process qq + T — > tt + T, and to repeat the reasoning of Jllf l by parametrizing the momenta of exchanged gluons 
along dijet total momentum.) The QCD factorization theorem |ll| predicts also that the interaction with the target 
via t-channel exchange by qq pair (which is expressed in terms of the skewed quark distribution) is not small at x > 0.1 
and moderate n t where gluon distribution is not large p6| , p6[ . At smaller x < 10~ 2 and n 2 > 1.5GeV 2 where the gluon 
distribution is large, as a result of x and Q 2 evolution, this term is a small correction to the exchange of the two-gluon 
ladder, see Eq. (|60|). So in this paper, we shall neglect this term. The two gluons are vector particles (bosons) in a 
color singlet state, so the dominant two-gluon exchange amplitude occurs in a channel which has positive charge and 
spatial parity, and is therefore even under crossing symmetry. Given this even amplitude, and the condition that we 
consider high energies v = 2p 7T mN and fixed small values of the momentum transfer t to the target, we may use the 
dispersion relation over invariant energy s at fixed t and fixed momenta of two jets to reconstruct full amplitude via 
the amplitude cut over intermediate states in s and u channels. In difference from the amplitudes of hard diffractive 
meson production, the amplitude for dijet production may have an imaginary part which varies with M 2 ^ et also. At 
the same time, within the approximations made here, the cut amplitude coincides with the imaginary part of the full 
amplitude. Furthermore, the relation discussed below makes it possible to reconstruct the real part of the amplitude 
from the imaginary part. For an amplitude corresponding to a slowly growing total cross section (A cx s a , a ~ 1) the 
relation is |27|| 

ReA(v,t) _ tt d lmA(u,t) 

W(i/,t) ~2dh^ n v ' (4) 

This means that we may simply calculate the imaginary parts of any contribution to the scattering amplitude, as a 
function of s and u, with the full amplitude obtainable from Eq.(|J). Furthermore, the imaginary part of the scattering 
amplitude, ImA varies rather slowly with is, leading to a small value of ReA/ImA. Thus ImA dominates in the sum of 
diagrams. The possibility of considering only the imaginary part of the scattering amplitude simplifies the calculations 
enormously. The relevant intermediate states are almost on the energy-shell (virtuality of quark, anti-quark, gluon 
<C k 2 ) and one can use conservation of four-momentum to relate the momentum of the relevant intermediate states 
to that of the initial state. 

There is another enormous simplification which is related to the issue of gauge invariance. The pion wave function 
is not gauge invariant, but the s,u cut parts of the amplitude 7r + g — > J J + g, for two gluons in a color singlet state, 
are calculable in terms of amplitudes of sub-processes where only one gluon is off mass shell. For such amplitudes the 
Ward identities J28|- the conservation of color current- have the same form as the conservation of electromagnetic 
current in QED. In QED the current conservation identity has long been used to simplify calculations of high energy 
reactions || . We will often use the Ward identities (2^] to extend the QED method to treat various contributions to 
our process. To be able to separate soft and hard scales one need to account for the cancellation of infrared divergences 
introduced by using the light-cone gauge A + = |^6| . Instead of accounting for this cancellation as in |4{J we choose 
different gauges for the description of parton distributions in a nucleon, and for the amplitude pion fragmentation. 
This is legitimate within the region of applicability of leading log approximations. Within the chosen gauge, the hard 
gluon exchange amplitude and fermion propagator have no infrared divergence. Note also that after demonstrating 
the factorization of the hard QCD amplitude from the soft QCD amplitude, we may and will approximate the soft 
part of the pion wave function by a system of free qq p7[ |. However, we stress that this approximation is dangerous for 
evaluating the pieces of amplitudes dominated by soft physics especially if propagators contain infrared divergences. 
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In that case, this approximation violates Ward identities and the energy-momentum conservation law. 

Additional common features arise from considering the relevant kinematics. In all of the two-gluon exchange 
diagrams we consider, Figs. 1-12, the target nucleon of momentum p emits a gluon of momentum k% and absorbs one 
of momentum k 2 ■ Conservation of four- momentum gives 

ki - k 2 =p-p' =Pf -Pit, (5) 

in which p' and pf are the final momenta of the target and the dijet. Taking the dot product of the above with p^, 
for the large pion beam momentum relevant here, leads to the relation: 

9 9 9 9 

mi — m f — 

x\-x 2 = —+ = — , (6) 



2pir -p 



where 



^1 2 

Xl, 2 = — , (7) 

and where rrif is the mass of the final 2jet system: m/ = M 2 j et . Within the parton model approximation 

x 2 > 0, (8) 

except the region of very small x 2 — denoted the wee parton region. In the parton model this condition follows from 
the requirement that a parton knocked out of a nucleon should be kinematically separated from the rest of the target. 
Otherwise the amplitude should be suppressed by a power of k 2 0. This suppression disappears at sufficiently high 
energies for which the parton wave function of a target develops wee partons. 

Another important consequence of the positivity of "mass" 2 of partons in intermediate states Eq. (|§[) is that the 
fraction j3 of the pion's (+) momentum carried by exchanged gluons should satisfy the condition: 

1 > > 0, (9) 

for our kinematics. The restriction (3 cx n"l/v can be justified within the leading lnl/a; approximation only. 

The results (|(||8|,||) are significant because they will be used in the evaluation of other diagrams. In particular the 
condition (|J) is not fulfilled for the diagrams in which the transverse momenta of quarks within the pion wave function 
are significantly smaller than the observed transverse momenta of jets. This means that the quarks in the pion must 
have very high transverse momentum to satisfy Eq. (|^). In this case, the quarks are closely separated and we may 
consider the configuration to be a point-like configuration. Such a restriction is operative in the kinematics for which 
the target wave function has no wee partons, i.e. for sufficiently large xn- On the contrary, for sufficiently small 
values of xjv , the sign of x 2 becomes unimportant because the amplitude does not depend on the sign of wee parton 
momentum. 

It is also worth emphasizing that the dominance of small size configurations in the projectile pion, so important 
to our analysis, is closely related to the renormalizable nature of QCD. This renormalizability implies, as extensively 
discussed below, that the selection of large transverse momentum final-state jets leads to a selection of the large 
transverse momenta of the quarks in the pion wave function, and also to some increase of transverse momenta of the 
exchanged gluons. 

One also needs to realize that the emissions in the in-state and absorptions in the out-state combine in calculating 
the usual parton density to produce the renormalized parton density. Thus it is necessary to guarantee suppression 
of gluon radiation collinear to the pion direction in the initial, intermediate, and final states. Otherwise an exclusive 
process will be additionally suppressed by powers of the Sudakov type form factor. This is a stringent condition which 
suppresses the contribution of all other diagrams except that of Fig. 1 because for small values of Xn, the time and 
longitudinal distance intervals (~ 1/ (2mjs[XN)) are easily long enough to accommodate the radiation of a gluon. If a 
pion is in a spatially small quark-gluon component, collinear radiation is suppressed because color is highly localized 
in the plane transverse to the pion momentum. As a result (similar to the case of meson production by longitudinally 
polarized photons) there is no Sudakov form factor type suppression for such processes [^9j . Note also that according 
to the QCD factorization theorem a pion in a small size configuration consists of qq pair accompanied by a coherent, 
relatively soft gluon field which follows the valence quarks without violation of coherence. This gluon field is included 
in the skewed gluon distribution. 

Consider now the impact of the above-mentioned condition for the interactions of a pion in a large size qq con- 
figuration. The q, q and gluons which start off far apart must end up with a qq pair close together in a final state 
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without collinearly moving gluons. In this case, the q and q must undergo a high momentum transfer without emitting 
gluons collinear to the pion direction. But such processes are well known to be exponentially suppressed by double 
logarithmic Sudakov-type form factors. Only in the case of a compact qq pair, of a transverse size commensurate 
with the virtuality of a gluon bremsstrahlung, would the gluon radiation be small. A related suppression, evaluated 
in Sect. Ill is the very small probability u>2 of finding a pion with q and q at average distances without a gluon field, 
if the probe has a resolution k\ 3> Aqcd- Note that under these conditions, in the typical parton configuration, 
gluons are experimentally observed to carry about ~ 1/2 of the pion momentum. Another example is pion scattering 
by a high momentum gluon field of a target. In the intermediate state there should be a strong collinear radiation 
along the pion direction because the color charge strongly changes its direction of motion, and there is no color charge 
nearby to compensate for this emission. This is similar to the effect of filling agap in the case of color unconnected 
hard processes like Higgs production via gg H in hadron-hadron collisions |3(J . 

In considering hard exclusive processes, one needs to address the problem of the end point contributions - the 
so called Feynman mechanism. We find that due to the color neutrality of the pion and the effect of target recoil, 
the amplitude for this mechanism is suppressed by a factor oc 1/k 2 as compared with that of the perturbative QCD 
mechanism. The contribution of the Feynman mechanism is also suppressed by powers of the Sudakov-type form 
factor and by the form factor W2- 

In previous papers |5|Jl7jl we have emphasized that the amplitude we computed in 1993 is calculable using pertur- 
bative QCD. However, there are five other types of contributions which occur at the same order of a s . The previous 
term in which the interaction with the target gluons follows the gluon-exchange in the pion wave function has been 
denoted by T\, However, the two gluons from the nuclear target can also be annihilated by the exchanged gluon 
(color current of the pion wave function). This group of amplitudes is denoted as T 2 . Another term in which the 
interaction with the target gluons occur before the gluon-exchange in the wave function of dijet has been denoted by 
T3. There are also terms, denoted as X4, in which the interaction with the target gluons spans the entire time between 
interactions with target gluons. This term corresponds to the interaction of the qqg configuration with target gluons. 
Still another amplitude, T5, describes the interaction of qq dipole with a target in non-leading order in a s ln«f / A-qcd- 

Here is an outline of the remainder of this paper. Sect. II considers the k\ dependence of the Feynman diagrams 
and selects Feynman diagrams having the minimal power of at fixed values of a s In k^/Aqcd in the lowest order 
in a s . We found that the processes where small size qq configuration is prepared before the scattering dominate 
diffractive dijet production. Subtle features of the arguments are discussed in Sect. Ill, which is concerned with the 
role of selection of exclusive processes in the suppressing of the contribution of hard processes related to inter-quark 
transverse distances 3> All of the terms T2-5 are shown to be negligible in the sense that they are smaller 

than Ti by at least a power of a s or by a factor of ^y, or by powers of Sudakov-type form factors and/or of a form 

factor W2, which is related to the probability of finding a normal-sized qq configuration of the pion. At the end of this 
section we demonstrate that the contribution to high n t dijet production resulting from the scattering of a large size 
qq dipole by a large transverse momentum (ss Kt) component of the target gluon field is suppressed at least by two 
powers of the Sudakov-type form factor and by the W2 form factor. Together sections II and III form the proof of 
the QCD factorization theorem for the diffractive dijet production. Rather general arguments for the small nature of 
the amplitudes T 2 _5 provided in the Sects. II and III are valid for the case of the photon projectile as well. This is 
because the contribution into the forward scattering amplitude of dijet production due to direct photon coupling to 
light quarks is oc bare mass of quark and therefore small |3l|]32[| . In the case of a charm quark photoproduction, the 
dominant term is given by a charm component of the direct photon wave function. 

Sect. IV is concerned with the evaluation of the dominant amplitude Ti, which has the form of the QCD factorization 
theorem in the leading order in a s and all orders in a s \n(K^r 2 ). The analysis performed in Sections II-IV shows 
that the z dependence of the leading (over powers of fcf ) term in the amplitude of diffractive dijet production is given 
basically by the factor z(l — z). The nuclear dependence of the amplitude, including a reassessment of the multiple- 
scattering correction of H, and nuclear shadowing effects is discussed in Sect. V. Experimental aspects, including 
the requirements for observing color transparency and the extraction of the coherent cross section, are discussed in 
Sect. VI. There is an electromagnetic background term, which becomes increasingly more important as n t increases, 
in which the exchange of a photon with the target is responsible for the diffractive dissociation of the pion. This 
process, which occurs on the nuclear periphery and is therefore automatically free of initial and final state interactions, 
is shown in Sect. VII to provide a correction of less than a few percent contribution to the cross section at values of 
Kt, s of the experiment |l4| ] but this correction rapidly decreases with an increase in the value of n t . A discussion of 
the implications of observing color transparency as well as a summary and assessment of the present work is provided 
in the final Sect. VIII. 
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II. SELECTION OF DOMINANT FEYNMAN DIAGRAMS FOR nN -» NJJ IN THE LEADING ORDER 

OF a s AND 

The kinematic constraints due to the energy-momentum conservation play an important role in the evaluation 
of amplitudes of diffractive processes. Therefore, we begin by deducing the necessary kinematical relations and 
introducing the light cone variables we use. Our interest is in the scattering at nearly forward angles. We denote 
momentum of the pion as p^, and that of the target as p. The three momentum of the nucleus in the final state is 
Pf z « M%j et /2Ei T = mpfX. The first relation is expressed in terms of the variables of the nuclear rest frame, and 
the second is in the variables of IMF of a nucleus where x — M^ et jv. The mass of the two jet system is given by 

M%j et — z (\_z) ~ K t/ Z (l — z )- We neglect the mass of quark m q as compared with the large jet momentum n t . The 
quantity v = 2(p 7r p)/A is the invariant energy of collision. Our notation is that z represents the fraction of the total 
longitudinal momentum of the beam pion carried by the quark in the final state, and 1 — z the fraction carried by the 
anti-quark. The transverse momenta are given by K t and —K t . A is the number of nucleons in the nuclear target. Our 
interest is in the kinematics for which the final state nucleus remains intact. This means that minimal momentum 
transferred to nuclear target —t m in ~ P 2 f z should be small: — t m i n R\/Z <C 1 i.e. x <C m ^ ~ A~ 1 / 3 /3. Here 
Ra — 1.1A 1 / 3 Fm is the nuclear radius. For small values of —t m i n R A /3, the effect of t m in can be easily accounted 
for because any form factor of the target can be approximated as exp tR\/6. 

For large enough values of n t , the result of the calculations can be represented in a form in which only the qq 
components of the initial pion and final state wave functions are relevant in Eq. ([l2]). See also the discussion at 
the beginning of section III. This is because we are considering a coherent nuclear process which leads to a final 
state consisting of a quark and anti-quark moving at high relative transverse momentum. It is necessary to examine 
the various momentum scales that appear in this problem. The dominant non-perturbative component of the pion 
wave function carries relative momenta (conjugate to the transverse separation between the q and q) of the order of 
Pt ~ , w 300 MeV. This is much, much smaller than the final state transverse relative momenta, which must 

be greater than about 2 GeV, the minimal value required to experimentally define a jet. The immediate implication is 
that the non-perturbative pion wave function, which is approximately a Gaussian, cannot supply the necessary high 
relative momenta. These momenta can only arise from the exchange of a hard gluon, and this can be treated using 
perturbative QCD. 

Restricting ourselves to Feynman diagrams having the leading power of s, «;£, a s (at fixed a s hiK, 2 / Aq CD and/or 
a s In 1/ x) and using a normal non-perturbative wave function which rapidly decreases with increase of the constituent 
transverse momentum gives the possibility of regrouping the diagrams into blocks having a rather direct physical 
meaning: 

M(N) = (T 1 +T 2 + T 3 + T 4 + T 5 ) . (10) 

Here the dominant terms T\ and T±b of Figs. 1 and 2 represent a type of impulse approximation, and we shall examine 
them first. Next we evaluate the possible role of color flow - of the interaction of target gluons with with gluon 
exchanged between quarks in the pion wave function -amplitude T 2 .This amplitude is also expressed through the 
same pion wave function as T\. The term T 3 corresponds to the final state interaction between jets. This will be 
followed by the discussion of the physical meaning of the other terms T±,T§ and the explanation of their smallness. 
At this point and below we rely heavily on the fact that, if a s -C 1 in the leading a s \nn 2 / Aq CD and/or a s hxl/x 
approximation, the sum of dominant diagrams (but not each particular diagram) has a ladder structure. In particular, 
various crossed diagrams which have a different form are needed to guarantee local gauge invariance for processes 
with large rapidity gaps, and to ensure the ladder structure of the sum of dominant diagrams. To derive such a ladder 
structure of the sum, it is important to explore crossing symmetry and the positive charge parity which follows from 
the dominance of the t-channel exchange two-gluon state of vacuum quantum numbers. Accurate exploring of gauge 
invariance is necessary for proving the ladder structure of the sum of Feynman diagrams. This structure makes it 
possible to include terms varying as lnl/a; and Iiik 2 /Aq CD , related to conventional QCD evolution, in skewed parton 
distributions or in skewed unintegrated gluon densities. In addition, we need to calculate the evolution of the pion 
wave function with transverse momentum and the interaction of gluons with this wave function. So we shall first 
classify and calculate skeleton diagrams, and then account for the QCD evolution. 

A. Dimensional estimate of the initial state hard interaction 

In our previous papers we investigated the term T\ of Fig. 1. QCD is a non-Abelian gauge theory in which 
an exchanged gluon may probe the flow of color within the pion wave function. So in the leading order of a s , 
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the requirements of gauge invariance mandates that we should consider also the related set of diagrams where an 
exchanged gluon is attached to the exchanged gluon in the qq component of the pion wave function, see Fig 2. 
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FIG. 1. A contribution to T\ a . The high momentum component of the pion interacts with the two-gluon field of the target. 
The displayed diagram occurs along with its version in which the gluons are crossed. Furthermore, there are four diagrams for 
each term because each of the gluons can be absorbed or emitted by either the quark or anti-quark of the beam pion. Thus 
only a single diagram of the eight that contribute is shown. 




FIG. 2. A contribution to Tib- The high momentum gluon current in the the pion wave function interacts with the two-gluon 
field of the target. The displayed diagram occurs along with its version in which the gluons are crossed, for different attachments 
of gluons in the pion wave function. Thus only a single diagram of the eight that contributes is shown. 



In the evaluation of the terms of Figs. 1 and 2, with color flow we use Gribov's observation Q that, within the leading 
a s In 1/x approximation, the polarization of gluons exchanged in the ladder is cx p where p is the four momentum of 
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the target. It is not difficult to show that the same gluon polarization dominates in the calculation of QCD evolution 
in the leading order in ln«;|/ Aq CD . It is also convenient to use the fact that only one gluon is off its mass shell (in the 
leading logarithmic approximation considered in the paper). This causes the equation that describes the conservation 
of color current to have the same form as in QED. We show here that the use of QED-type Ward identities ]2§| ] 
(allowed in computing the s, u cuts of the diagram, as explained in the Introduction) leads to the result that only the 
transverse components of the gluon momenta k\,k2 enter in the final result for the amplitudes Ti a ,Ti;,. Momentum 
factors, related to the contribution of the target gluons, are included by definition in the skewed gluon distribution. 

We examine the part of Ti a ,Tu> that arises from the exchange of the gluon k\ and that gives an on-shell qq 
intermediate state. The result is 

Tx = T la + T lb oc A^dTd^A^ (11) 

where A* ,N represents the gluon emission amplitude of the pion, nucleon and d^ v arising from the propagator of the 
gluons emitted or absorbed by the target. At high energies, the gluon propagator can be represented as M 



d »» K M _ (12) 
2ptt -p 



where pP is the nucleon momentum, so that 



Tla + Tlh0C 2 PK . P ■ (13) 

We will denote this presentation of the contribution of gluon exchanges-as the Gribov representation because he was 
the first to understand the dominance in the high energy processes of gluon polarization cx p, sec Q. Now we use 
current conservation, 

a; x • = o, (14) 

and employ Sudakov variables to describe the momentum k\ : 

kZ = x lp K + l3p»+k t . (15) 

We can determine the quantities a, f3 by taking the dot product of the above equation with either p or p^ and 
neglecting the relatively small factors of the square of the pion or nucleon mass. This gives 

k\ -Pir 

Xl = , 

P-Pll 

P=^. (16) 

P-Pir 

Using these results in the current conservation relation (Q) leads to the relation: 

A* ■» --La* •« A ^ x ' (17) 



The first and third terms of Eq. (|17| ) , in difference from the second, are proportional to s and therefore dominate 
over the second ||. Thus we find 

so that the exchange of the gluon k\ gives an amplitude (|T^ ) proportional to the small transverse momentum K\t- 
The net result of these considerations is that the contribution of the gauge invariant set of the diagrams including 
those of Figs 1, 2 takes the form: 

T la +T lb (X — ^- ^— ••• <y- r . (19) 

(2«i • Pit) (2p^ ■ p) 

The same trick can be made with the second exchanged gluon: 
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4A'\ • ki ut k 2 At 

r " + r - (a,^") -^-^- (20) 

The factors involving fc^ will be absorbed into the definition of the skewed gluon distribution of the nucleon. The 
useful result is in the denominator of Eqs. (p~9|j2C|) , because 



2ki • p = Xi s. (21) 



So 



*- +T »* (22) 

This formula in which the amplitude is expressed in terms of exchanges by transversely polarized gluons is the 
adjustment to QCD of the Weizsacker- Williams (WW) method of equivalent photons |33] ], of the Gribov derivation 
of reggeon calculus ||, and of the Cheng & Wu's impact parameter representation J34[ . So below we will denote 
such a formula as the Weizsacker- Williams representation of equivalent gluons. For the calculation of the dominant 
amplitude T\, such a trick is not very useful because effectively 12 < 11. So in this case we shall use the Gribov 
representation for the second exchanged gluon. 

Let us first perform power counting for the sum of terms T\ a + Tu within the Gribov representation for the 
amplitude. The number of strongly virtual propagators in Fig 1 is two and the number of large transverse momenta 
in the numerator from the vertexes is two. An additional factor 1/fcf follows from the cancellation of the sum of 
leading diagrams because of the color neutrality of the pion. In the following we include the factor a s //tf as part 
of the high momentum component of pion wave function. Slowly changing factors such as a sl lnx,hi(K 2 /Aq CD ) are 
included by definition in the skewed gluon distribution xG a and into the pion wave function. Finally we obtain: 
T\ a oc a s Xir(z, K$)xiGa{xii X2, K 2 )/k 2 . This result can be easily proved within the WW representation also. Note 
also that the comparison of the above result with the WW representation shows that the kinematical region around 
%2 ~ gives a negligible contribution to the integral over kit, P, x\. 

The contribution of the color exchange current - the term T^- is suppressed in the light-cone gauge as compared 
to Ti a by a power of Kt- To estimate the kt dependence of the diagrams, we shall use the WW representation of the 
gluon exchange between quarks in the initial pion and the target and the Gribov representation for gluon exchange 
with the target in the final state. Within this representation for the sum of diagrams, the cancellation between 

diagrams leading to the number of powers of \ is accounted for in a straightforward way. So the contribution of 

K t 

energy denominators gives: , L 3 . Additional oc A as compared to the product of gluon propagators is a result of 

cancellations between diagrams due to color neutrality of pion and final state of two jets. Account of gluon momenta 
in the nominator gives when both gluons in the gluon color current in the pion wave function are longitudinally 
polarized: (feit)rP/i,7rPA,jr$u a — 0. The contribution in the denominator when one of the gluons in the gluon color 
current in the pion wave function is longitudinally polarized, while the second is transversely polarized, also vanishes. 
This is because forward scattering cannot change helicity. The proof of this statement follows from the combination 
of the light cone gauge condition: = and the fact that the leading power of s is given by the + vertex for the 
interaction of aa nuclear gluon with a quark. Accounting for the color neutrality of the pion wave function and of the 
wave function of the dijct final state, as well as the anti-symmetry of 3 gluon vertex is also important. 

The use of identity: {kxt)r^ij,,t^\,t9 r u, \ — helps to prove that the leading contribution into Tu, is when both 
gluons in the gluon color current in the pion wave function are transversely polarized. Here k\ t is the transverse 
momentum of a gluon exchanged with target and K t is the transverse momentum of the jet. g r x is the Yang-Mills 
three-gluon vertex. 

The physical meaning of the obtained result for T\ a + Tu is that the interaction of quarks in the pion wave function 
with the target gluon with a relatively low virtuality is dominated by distances significantly larger than that involved 
in the pion wave function, and by the target gluon interaction with external lines of the amplitude without a gluon. 
This is the generalization to QCD of a theorem proved for QED by Low pq |. 

For completeness we present the result of the calculation of Feynman diagrams for T± in the leading a s In -np — 

QCD 

approximation where the integration over fraction of pion momentum [3 transfered to the target gluon by quarks in the 
pion wave function is not performed. The derivation is rather close to that in j36fl because in the leading a s In -r^ — 

H r QCD 

approximation the same polarization of target gluons dominates as in the leading a s In 1/x approximation. The 
deduced formula is actually very similar to the formula deduced in ]36| ] for hard diffractive vector meson production 
where in the leading In 1 jx approximation the integral over the unintegrated gluon density is replaced in the final 
step by the gluon density. So we will not repeat the detailed evaluation made in Q . 
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The result is: 



Ti/v = J a s (2(2 x (z, k 2 ) - x (z + (3, k 2 t ) - X (z - k 2 t )) + A( x (z + (3, k 2 t ) + X (z - (3, k 2 ))) ■ 

-F (3)7r — d kitjjrj' ( 2 ^) 



rf/K 2 ) 2 ^ /?" "k 2 k 2 2 
Here A is the two dimensional k t space Laplacian operator which acts on the pion wave function, fx = ~ ih ^^^w^ 



The first term of Eq. (|23|) is small and we shall explain why we neglect it. The factor: d ^ K a^ follows from the definition 
of the pion wave function and from the definition of hard amplitude in terms of series over the powers of the running 
coupling constant. The amplitude T± can be simplified by using the leading lnl/x approximation. In particular, one 
may express the amplitude in terms of the gluon distribution [p6| . This gluon distribution is however different from the 
usual gluon distribution determined from DIS processes because of the significant difference between the masses of the 
initial pion and final two jet systems. The necessary skewed gluon distribution is calculable for small values of x as the 
solution of QCD evolution equation, using the ordinary diagonal gluon distribution as the initial condition Here, 

in a fashion similar to [ p6| , we approximate: xG — J d 2 kit-grf5 fg /t, so that fx can be denoted the unintegrated 
gluon density. Equation ( |23| ) is a version of Eq. ([!]), in which the dipole approximation (keeping terms of the order 
r\ in the bracket [2 — exp(— i(rt ■ kit)) — exp(i(rt • kit))] in Eq. (M) is used. This approximation is reasonable even 
if (kitrt) is comparable or even larger than 1 (cf. discussion in |36| after Eq. (2.20)). Another useful approximation 
involves the value of f3 appearing in the argument of the pion wave function. The upper limit of integration over (3 
is dictated by energy conservation to be (3 < 1. For z ~ 1 one gets a further restriction that (3 < 1 — z. But in the 

2 

leading log -r^ 11 — approximation, the condition for the region of integration in (3 is more severe. It is given by the 

QCD 

requirement that the contribution of the target gluon longitudinal momentum into its propagator sa j3M 2 ^ et should 
satisfy conditions: (3M 2 j et w k\ t <C M| et . Hence in Eq.(23) it is legitimate to neglect (3 in the argument of the pion 



wave function, and to keep in Eq. fl23|) the term oc A only. The value of the lower limit of the integration: (3 > f3o is 
obtained from the energy conservation law, and from the QCD evolution which effectively suppresses the contribution 
of the region k\j (3qV ~ 1. 

Elastic interactions between high n t q and q in the final state may lead to an infrared contribution which is exactly 
canceled in the probability summed over the final state gluon radiation. It is important that the complete nature of 
final states allows the term T\ to account for the final state radiation and the space-time evolution of the qq pair in 
the final state. This phenomenon is a familiar feature of the theoretical analysis of the fragmentation of the small size 
wave packet, of e + e — > hadrons. 



B. Meson Color Flow Term-T 2 



The T2 or meson-color-flow term, of Fig. 3 arises from the qqg intermediate state, or from the attachment of both 
target gluons to the exchanged gluon or to the exchanged gluon and the quark(anti-quark). 
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FIG. 3. Contribution to Ti a of qqg intermediate state. The exchanged gluon interacts with with each of the target gluons. 
There is also a diagram in which the gluons from the target are crossed, and another two in which the exchanged gluon is 
emitted by the anti-quark. Only a single diagram of the four that contribute is shown. 

The diagrams of Fig. 3, considered in this section, have the same topology as the diagrams included in the term 

A 2 

Ti and therefore have the same power of Q f D , and lnl/x. In particular, they are included in T\ if l\ <C k%. (To 

K t 

visualize the relationship between T\ and T2 it is useful to move down the point of attachment of target gluons to 
the gluon in the pion wave function in the diagram Fig. 3.) In this case, a significant contribution may arise only 
from the perturbativc high momentum tail of the pion wave function. The effects of the gluon-gluon interaction may 
then be included as part of the target gluon distribution. In the case when the square of the transverse momenta of 

the gluon in the pion wave function is oc k% , this diagram can be considered as a non- leading order correction (to the 

2 

dipole-target interaction cross section) in a s at fixed values of a s In — . In contrast with the term T 2 , the definition 

QCD 

of T\ contains a requirement that the transverse momenta of gluons attached to quark lines should be much less than 
Kf, and that the pion wave function includes its perturbative tail. For the terms of Fig. 3, the non-perturbative pion 
wave function cuts off large quark momenta in the pion wave function. Thus gluonic transverse momenta are w k|, 
but transverse momenta of target gluons are still small: Kf t <C nf. The contributions of other diagrams with an almost 

A 2 

on shell qqg intermediate state, see Fig. 4, are not suppressed by a power of Q % D . 
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FIG. 4. Contribution to T2& from qqg intermediate state. The interaction of one target gluon field with exchanged gluon in 
the intermediate states. There is also a diagram in which the gluons from the target are crossed, and another group in which 
the exchanged gluon is emitted by the anti-quark. Only one of 16 diagrams that contribute is shown. 



It is easy to calculate the sum of the terms Ti a ^b in the framework of the Gribov representation where dominant 
contribution is given by the polarization of a target gluon which is oc target four momentum. Most straightforward 
calculation is for the ratio ^ because in this ratio all factors except the z dependence are canceled out. Really 
quark-gluon vertexes when gluons are attached to quark(anti-quark) lines in the pion wave function are the same for 
the terms 7\ , T 2 • In the light cone gauge all factors from the vertexes for the interaction of target gluons with a gluon 
in the initial and final state wave functions are effectively the same as for a gluon interaction with quarks except for 
the Casimir operator of the color group in the octet and the triplet representation. A subtle point of calculation is to 
evaluate the z dependence of this ratio. For certainty in the evaluation of term T 2 we assume that the nonperturbative 
pion wave function is equal to the asymptotic one. 

Thus the ratio is determined by the color content of color flow in the pion wave function and the quark color and 
by the dependence of energy denominators on the fraction of pion momentum carried by quarks and gluons. So 



T2 
Ti 



F 2 {8) 
f*(3) 



+ 



r(l-z) {l-zf 



\nz+^—^-\n{l-z)) 



(24) 



Here F 2 (i) (for i = 8,3 is the Casimir operator for octet and triplet representations of color group SU(3) C . The 
ratio H is w 0.5 for z = 1/2, remains nearly constant for \z — .5| < 0.3 and increases to 9/8 at z=0,l. This term is 
additionally suppressed by the Sudakov type form factor and by the form factor w 2 -see the discussion below. 



C. Final state interaction of qq Pair— T3 

The interaction with the target gluons may occur before the interaction between quarks in the final state, and the 
related amplitudes are denoted as T 3 , see Figs. 5 and 6. 
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FIG. 5. Contribution to T^ a - The high momentum component of the final qq pair interacts with the two-gluon field of the 
target. Only a single diagram of the eight that contribute is shown. 
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FIG. 6. Contribution to T^t- The gluon from the two-gluon field of the target interacts with the high momentum component 
of the final qq pair wave function . Only a single diagram of the eight that contribute is shown. 



The term T^ a includes the effect of the final state qq interaction. Fig. 6 includes the interaction of a target gluon 
with color flow in the wave function of final state. 

We need to evaluate only the s~,u— channel cuts of the diagram (and use Eq. ([|) to get any necessary real part). It 
is useful to define It as the quark transverse momentum within the pion wave function. Then there are two kinematic 
regimes to consider. The first has l t <§C K t , k lt <C n t , and the second l\ ~ k\ t ~ nf. We consider the former regime 
first, as it is expected to be more important. In this case, we shall employ conservation of the four-momentum to 
evaluate X2 ■ Conservation of the four-momentum can be used to relate the intermediate state (denoted by the vertical 
dashed line, occurring between the emission and the absorption of the gluons by the target in the diagram of Fig. 5 
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p7| ) of momentum p with p 2 = fn 2 with the intermediate state. The mass of the qq intermediate state is given by 

fh 2 x\v — x\(5v — k 2 t , (25) 

— k~ — k~ 

where (3 is the light-cone fraction of the pion momentum carried by an exchanged gluon: (3 = = Thus we 

arrive at the equation: 

It follows from the requirement of positivity of energies of all produced particles in the intermediate states that 
< j3 < 1. We can now calculate m 2 directly in terms of the light cone momenta of the qq pair in the intermediate 
state: 

m 2 =(| + J^Z^)(i- / 3)_ j £. (27) 



z 1 — p — z 



Combining Eqs. (E5j),(E7n we obtain 



which, when using Eq. ^) leads to 



„ + filial!) = ^+*„ <2 9 ) 

v V z 1 — p — z 1 v 



Therefore 

z 1 — — z z(l — z) 

In order for the term T^ a to compete with T± a we need to have l f < k ( , k\ t <C n t -otherwise T^ a will be additionally 
suppressed by the power of n 2 ,a s . These kinematics cause Eq. ( |30| ) to yield the result: — X2 <x k 2 /v. 

This argument can be carried out for all combinations of diagrams represented by Fig. 5. For example, another 
attachment of gluons, in which the gluon k\ is absorbed by the quark, corresponds to interchanging z with 1 — z, 
and therefore leads to the same result for x-i- Evidently this result for X2 is valid in the leading a s \iiK 2 /AQ CD 
approximation also. Thus we consider the second situation: If ~ k\ t <~ k 2 . In this case, the initial pion wave function 
contains a hard quark, and we discuss hard radiative correction in the next order of a s . This is the typical situation in 
which there are extra hard lines, as compared with the dominant terms, and one obtains a suppression factor ~ 
which could be compensated by d 2 kt integral. However this integral does not produce Itik 2 /Aq CD because the region 
of integration is too narrow. So this contribution is at most NLO correction over a s . But we restrict ourselves by the 
LO contribution only. 

The presence of color flow in the wave function of the final state leads to the interaction of a target gluon with a 
gluon in the wave function of the final state, see Fig. 6. This term is suppressed by an additional power of 1/k 2 - The 
proof of this statement repeats the same reasoning as that explaining the suppression of the term . It heavily 
uses the WW and the Gribov representations, discussed in subsection A, and the identities which follow from the 
antisymmetry of the vertex for the three gluon interaction, the color neutrality of the pion wave function, and the 
dijet final state. In the derivation it is helpful to use the observation that effectively |a^2 1 oc >hI v - Evidently similar 
reasoning is applicable in computing amplitudes to leading order in a s and all orders in a s \n.n 2 /l 2 . 

Repeating the same reasoning as in the estimate of the terms of Ti a ,Tu,, and remembering that —X2 ~ 4l v we 
achieve the estimate: T3 oc a 2 xiGA(xi, x%, K 2 )/K,f. It is instructive to investigate whether the Feynman mechanism, 
where the leading quark (anti-quark) carries a fraction of the pion momentum z' close to 1 but high momentum jets 
are formed by the action of a final state interaction, may compete with the PQCD description. In this case transverse 
momenta of constituents l t in the pion wave function are expected to be equal to the mean transverse momenta of 
partons in the non-perturbative regime. For certainty let us model the Feynman mechanism by assuming that recoil 
system is quark(anti-quark) with momentum \ — z' close to . Within this model we will obtain Feynman diagrams for 
the term T2, but with the region of integration defined by the Feynman mechanism. A simple dimensional evaluation 
of term T3 due to the Feynman mechanism within the Gribov representation shows that it is suppressed by the 
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powers of n t . The contribution of the region ppr^) "C -^fjet nas been considered above-it is additionally suppressed 
for Feynman mechanism by the restriction of the region of integration over z' . Thus our next discussion is restricted 
by the consideration of the contribution of the region: (j-fjry > Affjet . 

T >^J*' l ^M^,<B) d%d " <31) 

In the above formulae we use Brodsky-Lepage convention for the definition of wave functions and retain terms 

i 2 

maximally singular when z — > 1. Power counting is simple: the factor /-,_'/■, 2 is from the gluon exchange in the 

final state. The factor I 2 /(I — z 1 ) is singular when z 1 — * 1. It originates from the quark vertexes accompanying 
the propagator of the gluon exchanged in the wave function of final state. Here 1/(1 — z') follows from a transition 
when a fraction of the pion momentum carried by quark tends to 0. The factor M 2 (2jet) — l 2 /z'{l — z') in the 
denominator is due the fermion propagator adjacent to the hard gluon exchange in the wave function of the final 

state. Here Mf nt « m [" z i * is the mass of an intermediate state, and m 2 ec is the invariant mass of the recoil system 
in the Feynman mechanism. In the region of integration 1 — z' -C nf-^ijet one ma y neglect by M 2 (2jet) in the 
denominator as compared to l 2 /(l — z'). So one obtains : T2 oc J ^(z' ,l 2 )d 2 ltdz' . In this case, another factor 

of 1/kJ arises from the integration over z'. Hence we have found that the Feynman mechanism is a higher twist 
correction to the PQCD contribution. The Feynman mechanism is further suppressed by the requirement of a lack of 
collinear to pion momentum radiation-see the discussion below. 



D. Gluon admixture to the wave functions of initial and final states - T4 



The Feynman diagram corresponding to Fig. 7 contains the time ordering corresponding to the qqg configuration in 
the pion wave function interacting with the quarks in the final state. In taking the imaginary part of the amplitude, 
the intermediate state must contain a hard on-shell quark and a hard on-shell gluon. But such a state cannot be 
produced by a soft almost on-shell quark in the initial state, so there is an additional suppression factor, caused by 
the rapid decrease of the non-perturbative pion wave function with increasing quark virtuality. This factor is greater 
than a power of k 2 . One may also consider the case when the transverse momenta of quarks in the pion wave function 
are large enough to use PQCD. Then the large virtuality of the quark introduces a suppression factor of -^7, with at 

least one power of arising from the quark line for the transition q — > qg and another factor of A arising from the 

K t 1 L t 

pion wave function. There are additional factors: \ arises from the hard fermion line, and from the application 

K t K t 

of Ward identities and the condition: x\v, x^v oc n 2 . A factor of k 2 1 2 is present in the numerator, with k 2 originating 
from the vertices in the WW representation and l\ from the integration over quark momenta in the pion wave function. 
All in all this amplitude is suppressed by the factor l 2 /(n 2 ) 3 . Another case occurs when I 2 oc n 2 . Then this diagram 
will be suppressed as compared to T\ at least by one power of a s without the large factor ln«;|/Ag C£) . But here we 
restrict ourselves to the analysis of LO corrections. 
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FIG. 7. A time ordering that contributes to T4. The qqg state interacts with the target. Only a single diagram of the eight 
where a gluon interacts with quarks in a pion fragmentation region that contribute is shown. 

Similar reasoning helps to prove that the contribution of diagrams in Fig. 6 is suppressed by a factor l\j as 
compared to that in Fig. 1. This is the power-type suppression if the pion wave function is non-perturbative, and 
may be a NLO a s correction if the perturbative high momentum tail is included in the pion wave function. 




I 

FIG. 8. A contribution to T^. The target gluon absorbs a gluon of pion wave function. Only one diagram of the eight that 
occur is shown. 

Another contribution to T4 arises from the sum of Feynman diagrams in which the gluon exchange between the q 
and q in the beam occurs during the interaction with the target, see Figs. 8, 9, and 10. The naive expectation is that 
such terms, which amount to having a gluon exchanged during the very short interaction time characteristic of the 
two gluon exchange process occurring at high energies, must be very small indeed. 

The intent of this sub-section is to use the analytic properties of the scattering amplitude to show that T4 is 
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negligible. Instead of calculating the sum of the imaginary parts of all of the amplitudes, we will prove that this sum 
vanishes by analyzing analytic properties of the important diagrams. Each considered diagram contains a product 
of an intermediate-state quark and anti-quark propagator. At high energies, these propagators are controlled by the 
terms of highest power of x\ 2p-p 1T — xiv, and (as to be shown) have poles in the complex x\ plane which are located 
on one side of the contour of integration. The sign of the term containing (u) in each propagator unambiguously 
follows from the directions of pion and target momenta. If we can show that the typical integral is of the form 



dxx- — ±- — -, a,/3>0 (32) 

[a x\v — a + ie)\p X\v — b + ie) 



the proof would be complete. 




FIG. 9. A contribution to T± c . The quark absorbs a gluon of momentum fci, exchanges a gluon with the anti-quark, and 
then emits a gluon of momentum fe. Only one diagram of the eight that occur is shown. 



We now consider the Feynman graphs, starting with Fig. 9. Once again we compute the imaginary part of the 
graph and consider the intermediate state as being on the energy shell. The propagator for the line (a) has the factor 

(fci + z'pw) 2 - ml — z'xiv + • • • , (33) 

while that of the near-mass-shell line (b) is independent of Xi, because the quark momenta in the final state and in 
the pion wave function are not connected with the target momentum. The propagator of line (c) has the factor 

(fe + <7i) 2 — m 2 q = X2 z v + ■ ■ ■ = x\ z v + ■ ■ ■ . (34) 

Here q\ is the momentum of thejet (z, Kt) and • ■ • denotes the terms which are independent of x\. The last equation 
is obtained from using Eqs. (§,|). The results (H||) show that the dia gram of Fig. 9 takes on the mathematical 
form of the integral (pi). Thus this term vanishes. 
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FIG. 10. Another contribution to T4 — T^a- The quark absorbs a gluon of momentum fci, exchanges a gluon with the 
anti-quark, and the anti-quark emits a gluon of momentum fe. Only one diagram of the four that contribute is shown. 



We also consider the diagram of Fig. 10. In this case there are three propagators (a), (b), (c) that have a term 
proportional to X\V, but the coefficients are not all positive. The propagator factor for line (a) is given by 

(x% ^ + fci) 2 = x\z'v + ■ ■ ■ , (35) 

while that of line (c) is given by 

(fc 2 + q 2 ) 2 — to 2 , = £2(1 — z)v + • • • = (1 — z)X\V + ■ ■ ■ (36) 

At the same time, the coefficient multiplying xi in the propagator (b) (gluon production) has no definite sign. Thus 
for this diagram the integral over x\ does not vanish. The presence of an additional gluon in the intermediate state 
means that x\ cx x%. This and the use of the Weizsacker- Williams representation allow us to conclude that this 
diagram is suppressed by a power of k 2 as compared to T\. Similar logic can be applied to any of the diagrams 
contributing to T4. The physical idea that the intermediate qq state does not live long enough to exchange a gluon is 
realized in the ability to close the contour of integration in the upper half plane or in the suppression by the power of 
Kq . The analyticity of the scattering amplitude in the upper half is a consequence of causality. Thus the physical and 
mathematical ideas behind the vanishing of T4 are basically equivalent for all diagrams at high enough beam energies. 

We conclude this section with a brief summary. We analyzed the leading diagrams for the pion dissociation into 
two jets and found that the hard dynamics-amplitude T\ determines this process, with the initial state wave function 
determined by the hard gluon exchange diagram. In the next section we shall show that amplitudes T%, T3 are strongly 
suppressed by the requirement of the lack of radiation in the final state. This is because in the lowest order in a s these 
amplitudes correspond to the propagation of a nonperturbative qq dipole or a q q g tripole over large longitudinal 
distances. 

While we were revising the manuscript in response to referee's questions, a paper |38| appeared which claimed that, 
if one includes terms beyond the leading-logarithm approximation in a s lnx, factorization does not hold and end-point 
singularities break collinear factorization. Our calculation shows that such problems are not present in the leading 

fe 2 

order approximation which keeps terms leading order in a s In — . The use of our light cone gauge, where the 

QCD 

asymptotics of the fermion propagator has no infrared singularities in the hard regime, makes the separation of scales- 
QCD factorization- rather straightforward. Moreover in this gauge, terms ~ lnl/x in the hard regime are related 
to the exchange by the gluon in the multi-Regge kinematics only Q. As a consequence of the QCD factorization 
theorem terms ~ lnl/a; are included in structure functions of the target. On the contrary, if one uses the standard 
gauge A + = or Feynman gauge, the cancellation of infrared singularities occurs if one includes the renormalization 
factor arising from the hard Fermion propagator, cf. discussion in Ref. [46|. Taking these terms together, the resulting 
amplitude contains no end-point singularity, and the factorization holds. 
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III. POST-SELECTION OF THE PROJECTILE WAVE FUNCTION BY FIXING THE FINAL STATE 



A specific feature of the processes considered here is that the restriction on the composition of the final state selects 
a rather specific initial configuration of the projectile hadron. Following S.Nussinov (to be published) we denote such 
measurements as post-selection. In our case, the initial and final state wave functions are built at large longitudinal 
distances oc 1/2itinx, so there is plenty of time for radiation to occur. But our definition of the final state forbids 
radiation collinear to the direction of pion momentum. In particular, the relation between the transverse momentum 
of jet Kt and the mass of the diffractively produced system, M 2 (2jet) = — z), is natural only for a qq final 

state. Thus any contribution due to processes for which such radiation is kinematically and dynamically permitted 
must be suppressed by the powers of Sudakov-type form factors and by a W2 form factor. The first step is to analyze 
how the trigger for the two jet state restricts the composition of the final state. 

A. Three jet production. 

A question arises^ whether the trigger used in Ref. Jl4| allows the separation of production by a qq state from 
the production by a qqg state (and more complicated states containing relatively soft partons) as the source of the 
observed dijets, and whether the presence of such states may change the n t dependence of cross section. Such states 
are certainly important in inclusive diffraction in DIS at HERA. A typical diagram corresponding to such a process 
is presented in Figure 11. 

9 ^ ^ Z(l-Z 2 ) 
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FIG. 11. A typical diagram corresponding to production of large mass qqg state. 

(l — z )M 2 - -\-k 2 

An analysis of kinematics shows that the virtuality of the gluon interacting with the target is « - — — — — 

where z g is the fraction of the pion momentum carried by a gluon in the final state and k tg is its transverse momentum, 



and as usual, the mass of the final dijet system is given by = M%j et = k^/z(1 — z). The mass of a qqg system in 

the final state is given by M$ jet = ff^ ^ + k^ g /z g . If z g < k^ g /M^ jet then M| je4 > 2M$ jet , and therefore the qqg 
state is distinguishable from the two-jet state by the relation between k\ and the total mass. In this kinematics, the 
cross section depends on Kt as a s 

( ffi3 (".(^>ffii.*».*?)) 3 j Beyond this kinema tics, the cross section is additionally 



1 We thank J.Bjorken, and D.Soper for asking this question. 
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suppressed by the factor cx {z g M 2 ,^ et )~ 2 . Thus a contribution of such configurations is either dynamically suppressed, 
or distinguishable from the contribution of a qq configuration. Note also that the feature that k tg is small would not 
lead to the cancellations needed for color transparency phenomenon to occur, and therefore the A-dependence would 
differ from what is observed. 

It is important to be able to distinguish dijet production from the soft diffraction into a hadronic state having a 
total mass and containing leading twist dijet production with mass M| jet . Evidently Mf^ < M^f*. Cross 

section of such diffractive processes are often described as the hard diffractive mechanism of Ingelman-Schlein [fj9| . 
In this mechanism one considers a hard scattering of a parton belonging to the projectile it with a light-cone fraction 
x„ off a parton belonging to a diffractive parton density (effective Pomeron) with a light-cone fraction (3xp . 

Here 1 — xp is the nucleon momentum in the final state, cf fl39f . It follows from the above definitions and energy- 

m 2 

momentum conservation law that i?" 1 = . A 3 ' 1 = x„f3 Thus Rr 1 is the fraction of the total mass of the diffractive 

M <uff 

state carried by two jets. We are interested in the paper in the limit when the mass of a diffractively produced 
hadronic state is carried mostly by 2 jets i.e. when R = Rq ~ 1. In this case the cross section for the production 
of dijets with R > R (where R being close to 1 is determined by the accuracy of measuring a fraction of the pion 
momentum carried by two jets) is 




Here a factor k 4 is the usual k dependence of the cross section of the hard two parton collision, yjt — is the usual 

M diff 

Pomeron flux factor, the factor (1 — x^) 2 is the parametrization of the parton density in the pion at x„ — ► 1, fp(/3) 
is the diffractive structure function. Taking fp{(3) oc (1 — (3) for /3 — > 1 (see e.g. pQ] ) we obtain da/dn 2 {R > Rq) cx 
(1 — Rq) 5 /k 6 . The factor (1 — i?o) 3 arises from the integration over x^ and factor (1 — Rq) 2 is from the integration 
over f3. The additional factor is because in the discussed region s» Mfj et . There is also a contribution of the 

Coherent Pomeron [[0J corresponding to j3 = 1. It leads to a similar suppression as a function of Rq: cx (1 — i?o) 4 - 

Hence we conclude that the importance of the discussed leading twist mechanism as compared to the the exclusive 
dijet production term which is cx k~ s depends on the degree of exclusiveness of the experiment. The experiment 
]14j imposes the condition that M^jet = Mdiff- Such a selection should have been rather efficient since due to the 
acceptance of the detector a condition was imposed that all produced pions should have momenta larger than a 
minimal one. 



B. Suppression of the final state interaction 

Let us consider the important consequences of the formulated above restrictions on the composition of the final 
state. The existence of the term T 3 displayed in Fig. 5, caused Jennings & Miller |42| to worry that the value of Mn 
might be severely reduced due to a nearly complete cancellation. However, we shall explain here that this term as well 
as the term T2 are strongly suppressed in QCD as compared to the naive PQCD calculation explained in Sect. II. This 
suppression is the only way to resolve an evident contradiction: the kinematic restriction on the final state discussed 
above forbids radiation collinear to the pion momentum, but such radiation naturally arises in a hard collision, or from 
the presence of a significant gluon admixture in the non-perturbative pion wave function when nonperturbative qq and 
more complicated configurations propagate large and increasing with energy longitudinal distances. It is convenient 
to represent this suppression as a product of two factors: w — w\W2- The first factor accounts for the well understood 
suppression of the collinear initial state radiation in the scattering of a target gluon off a low k t quark. Remember 
that by definition a non-perturbative pion wave function does not include a PQCD high momentum tail. Accounting 

k 2 

for the LO QCD evolution will not change this conclusion in the LO approximation over parameter a s In -rj- 4 — . In 

A QCO 

the case of the amplitude T3 the diagrams where the hard gluon exchange is present both in the initial and the final 
states are potentially dangerous. However, this contribution into T 3 is suppressed by the power of a s as compared to 
the amplitude T\ . Similar radiation is permitted for the processes described by the amplitude T3 because tripole qgg 
propagates large longitudinal distances. This radiation in the direction of pion momentum carries a finite fraction of 
pion momentum which is significantly larger than that for wee hadrons which are products of the jet fragmentation. 
These wee hadrons carry cx {fpj^ fraction of pion momentum. The w\ form factor suppression for the amplitude T3 

is given by the square of Sudakov type form factor: S 2 (k 2 /If); one form factor arises for each collision with a target 
gluon. Here l t is the transverse momentum of a quark within the non-perturbative pion wave function. This form 
factor is a square root of the form factor of quark ||3| because the radiation of gluons off the final quark is included 
in the definition of the final state. In the light cone gauge: 
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w x = S 2 {k 2 /1 2 ) = exp(-4/3gln 2 K 2 /* 2 ), (38) 

where, using the leading log approximation, we replaced k\ t by Kt in the argument of the Sudakov form factor. Recall 
that the form factor S 2 is exactly the Sudakov form factor which enters in the Dokshitzcr-Gribov-Lipatov-Altareli- 
Parisi evolution equations as a coefficient of the S(x — 1) term, see the discussion in Ref. Q. 

w\ form factor should be practically the same for both amplitudes T2 and T3. This is because of total transverse 
momentum of the system of high k t quark(antiquark) + high k t gluon in the amplitude T3 is small and controlled by 
the nonperturbative pion wave function. Therefore in the collision with a target gluon this quark(antiquark) - gluon 
system radiates in the direction of pion momentum as a single quark(antiquark). 

The second suppression factor, W2, accounts for the dependence of the form of the effective QCD Lagrangian and 
appropriate degrees of freedom on the resolution. As a result of partial conservation of axial current f v is independent 
on the resolution-invariant on renormalization group transformations. At the same time important degrees of freedom 
depend on the resolution: dominant degrees of freedom in hard processes are bare quarks. On the contrary, in the 
non-perturbative regime because of the effects of spontaneously broken chiral symmetry dominant degrees of freedom 
are constituent quarks, pseudogoldstone modes -pions, various condensates of quark, gluon fields, etc. Remember 
that for the average-sized configuration of a hadron, approximately half of the pion momenta is carried by gluons. 
(Within a constituent quark model renormalizability of QCD is usually accounted for by introducing effective mass and 
effective interaction for the constituent quarks and is calculated in these models in terms of the constituent but not 
bare quarks. The evaluation of /„- within a bag model should include a prescription how to treat the nonperturbative 
volume energy density - bag surface. Within the Weizecker- Williams approximation this energy is equivalent to the 
gluon cloud in the light cone wave function of a pion .) 

In the case of interaction of local current f n is calculable in terms of the distribution of bare quarks which accounts 
for both nonperturbative and perturbative effects . In this case the suppression factor u>2 tends to one for A; 2 , — > 00 
since in this case only short distance perturbative degrees of freedom survive. So amplitude T\ containes no additional 
suppression factor W2- On the contrary in the amplitudes T2,T% the structure of constituents in the non-perturbative 
pion distribution is resolved as a result of a large time interval between two consequent hard collisions of pion 
constituents off target gluons. A hard collision of the pion's constituents with target gluons (with momentum transfer 
(k 2 <C k 2 )) frees some gluons from the pion wave function, and emitted radiation is collinear to the pion momentum. 
This radiation is, however, forbidden by the restriction on the final states discussed above. Thus the selection of a 
component of the pion wave function of a size determined by non-perturbative QCD, leads to an additional factor W2 
which suppresses this contribution to the scattering amplitude. 

We may estimate the factor W2 using models as follows: 

= sj dj* j_i ggft-g - z 2 re(6-i +Z1 + Z2M1- z 2 - go 

f^dztdzzil- z t - z 2 ) n 6(l- Z!- z 2 ) v a; v j 

Here are the fractions of the pion momentum carried by jets in the final state, and S is the experimental 
uncertainty in the fraction of pion momentum carried by high re t dijets. Within the democratic chain approximation, 
which reasonably describes the Q 2 and x dependencies of hard exclusive processes, the value of n = 1 for the 
component which contains qq and a valence gluon. Another estimate of n can be obtained in the constituent quark 
model assuming, for simplicity, that each constituent quark consists of a bare quark and one gluon. In this case 
we are effectively dealing with a four particle system and hence n > 2. The factor W2 gives a significant numerical 
suppression for all the models. No such suppression appears in the amplitude T± because the radiation of a gluon with 
small transverse momentum from a highly localized (e.g. size l/ft*) color-neutral quark configuration is suppressed 
by a power of k 2 . (Note that in the case of the final state interaction, the non-perturbative pion wave enters at 
inter-quark transverse distances. This means that gluon radiation with transverse momenta, k 2 g < I 2 , should be 
suppressed. For gluon radiation with larger transverse momenta there is no restriction.) Amplitudes T2, T3 correspond 
to the propagation at large longitudinal distances of "large" size qq pair whose transverse size is controlled by non- 
perturbative pion wave function The net result of all of this is that amplitudes T27 T3 and Figs. 3,5,6 can be neglected 
and we shall therefore ignore the amplitudes T2 , T3 . 

We also observe that the contribution of Feynman diagrams in which the skewed distribution is modeled by the 
scattering of an on-shell quark or gluon is suppressed even further by the square of the Sudakov form factor, in 
addition to form factors Wu w 2 . 
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C. Non-leading order approximation to the dipole-target interaction 

In the previous section we used the leading a s In re 2 /A-q CD approximation (for the case k\ 2t ^ K t) to estimate 
diagrams. To go beyond this, one should take into account the contribution of the region kf t as k\ to obtain the 
qq dipole-target interaction. The related contribution to dijet production is denoted as Tg. To some extent this 
contribution has been discussed above, in relation to the amplitude Tj. For the amplitude Tg, quarks in the wave 
function of the initial pion have small transverse momenta /(, and there is no hard interaction between the q and q in 
the initial or final states, see Fig. 12. 
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FIG. 12. T5, Hard color flaw diagram. 

Let us outline various phenomena relevant for the smallness of this contribution as compared to T\. 

1. The amplitudes given by diagram Fig. 12 correspond to two hard collisions occurring at the different space-time 
points where color abruptly changes the direction of its motion. The longitudinal distances (time interval) 
between both hard collisions are large, as 2Et T /k%, so that the emission of gluon radiation is permitted. The 
difference in impact parameters characterizing both hard collisions is as , where l t is the transverse momentum 
of quark in the non-perturbative pion wave function. Thus collinear gluon radiation in both collisions is not 
canceled, and forward gluon bremsstrahlung (precluded in T\ by the localization of color in transverse space) 
occurs. Such radiation effects are precluded by our choice of final states for which radiation collinear to pion 
momentum is absent. Thus the ratio Tg/Ti is suppressed (as explained in the previous subsection) by the 
square of a Sudakov-type form factor: S 2 (k^/1^) — exp(— 4/3|^ In 2 /e 2 t /Z 2 ). Including the effects of the Sudakov 
type form factor leads to an increase of the effective value It with an increase of k\t, and therefore will change 
kinematics of this diagram in the direction to that for the term T\ . 

2. Final states with additional gluons collinear to pion momentum are initiated by the components of the pion 
wave function which are close to the average. The significant probability for such processes follows from the 
fact that the gluons carry ~ 1/2 of the pion momentum. Thus the selection of the component of a pion wave 
function of average size, but having no gluons, leads again to the appearance of the additional suppression factor 
u> 2 . 

3. As the consequence of the rapid decrease of the non-perturbative pion wave function with l\ for the diagrams of 
Fig. 12, the target gluon actually has a negative value of %<x. To see this, we apply Eq. ([30j) to the situation when 
the transverse momenta of quarks in the initial state It are much smaller than Kt, so that kit ~ Kt, kit ~ —Kt- 
In this case 
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Here (3 is the fraction of the pion momentum carried by a target gluon. To satisfy the condition X2 > one 
needs 

1 - z > (3 > (1 - z) 2 , or z > (3 > z 2 . (41) 

This condition cannot be satisfied within the leading logl/x approximation where (3-^1. Moreover, Eq. ( |40| ) 
implies that — x-iv oc k 2 , except for a narrow interval in (3 near (3 = (1 — z) 2 . The contribution of this 
narrow interval is suppressed by the small length of this interval. (The end point z ~ 0, z ~ 1 contributions 
are suppressed by the pion wave function.) In the non-perturbative regime- within the parton model- the 
contribution of the Xi < region cannot be expressed in terms of the unintegrated gluon density, and it is 
suppressed by power of 1/n 2 Q|. On the contrary, within the leading logl/x approximation such a strong 
suppression is substituted by a slower one whose value depends on the ratio of k 2 to the square of the mass 
of the recoil system (the wave function of a target where a high K t gluon is removed). This ratio is k 2 /(3v. 
Evidently, the contribution of small j3 is suppressed by the rapid decrease of structure function with k 2 . This 
cutoff is significantly stronger than that for the term T\ because of smaller transverse momenta of target gluons. 

Thus we conclude that the selection of the final state leading to the post-selection of the initial state shows that 
approximation of having a point-like pionic configuration in the final state should be valid even beyond the leading 
order approximation used in this paper. Therefore formulas representing a hard process as the convolution of a pion 
wave function with the interaction cross section || should be valid beyond the LO approximation. 

It has been suggested in [[IB) that gluon scattering off on-mass shell quarks and gluons in the target wave function, 
see Fig. 12, is the dominant process for the diffractive dijet production by a nuclear target. We explained above that 
in QCD this amplitude is suppressed as compared to the term T\ by the product of the square of Sudakov type form 
factors and the form factor wi. So we will neglect now and forever the amplitudes presented in Fig. 12 . 



IV. AMPLITUDE FOR nN -> N JJ-EVALUATION OF THE DOMINANT TERM 

Let us consider the forward (t = t rn i„ sa 0) amplitude, M., for coherent dijet production on a nucleon ttN —> NJJ 
§: 

M(N) = (f(K U z),N'\ ?\ir,N), (42) 

where / represents the interaction with the target nucleon. The initial \ir) and final \f(nt,z)) states represent the 
physical states which generally involve all manner of multi-quark and gluon components. For large enough values 
of n t , the result of calculations can be represented in a form in which only the qq components of the initial pion 
are relevant in Eq. (f42|). This is because we are considering a coherent nuclear process which leads to a final state 
consisting of a quark and anti-quark moving at high relative transverse momentum. 

We showed above that the dominant (in powers of s and n 2 ) Feynman diagrams will be expressed in terms of the 
light-cone wave function of the pion. We therefore begin by calculating the high transverse momentum component of 
the pion wave function. The non-perturbative component of the light-cone pion wave function is represented by | tt) , 
and the high momentum components can be treated as arising from the following approximate equation p6j : 

I n qq ) = G V eff | tt), (43) 

where Go(pt,y) is the non-interacting qq Green's function for p 2 3> m 2 ,m 2 r 

(p t ,y\G \ P W) = Si2)(Pt - p ^ {y - y '\ (44) 

2/(i-y) 

in which m q represents the quark mass, y and y' represent the fraction of the longitudinal momentum carried by 
the quark; and the relative transverse momentum between the quark and anti-quark is pt- The complete effective 
interaction, obtainable in principle from PQCD, implicitly includes the effects of all Fock-space configurations. 

The evaluation of the graphs corresponding to Fig. 1 consists of two parts. As a first step, the Feynman diagram 
of Fig. 1 can be rewritten as a product of a high momentum component of a light cone pion wave function with the 
amplitude for the scattering of a quark-anti-quark dipole by a target. So, we need to know the relevant part of the 
pion wave function. Secondly, we need to determine the interaction with the target (here with the gluon field of the 
target) which causes the pion to dissociate into a qq. 
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A. High momentum component of the light cone pion wave function^ 

The full wave function | tt) is dominated by components in which the separation between the constituents is of the 
order of the diameter of the physical pion. But there is a pcrturbative tail in momentum space which accounts for the 
short distance part of the pion wave function. This tail is of dominant importance here because we take the overlap 
with a final state constructed from constituents moving at high relative momentum. It is therefore reasonable to start 
our considerations from the one gluon exchange contribution V 9 , and use the light cone gauge: p' v ^A^ = and the 
Brodsky-Lepage | j46[ normalization and phase-space conventions. Here p'^ = p-^ — c-pw where c is determined from the 
condition: p'^ 2 = 0. Evidently for sufficiently high energies of pion projectile such a gauge will almost coincide with 
the gauge A_ = and therefore the Wilson line operator between qq pair: exp(i J A^dx^) ~ exp(i J A t dx t ). Thus in 
the target rest frame, at the light cone where X- = the Wilson line operator does not produce additional gluons in 
spite of the large and increasing with energy longitudinal distances x+. The wave function x{kt, x) is gauge invariant 
because it depends on transverse, i.e. on physical degrees of freedom. (On the contrary in the Feynman gauge the 
qq pair evolves into a many particle state because of the Wilson line operator. So a calculation in the Feynman 
gauge should include the evaluation of the form factor which guarantees dominance of the two particle final state.) 
Another advantage of the gauge chosen in the paper is that the amplitude for hard gluon exchange has no infrared 
divergences and therefore the separation between large and small distances is straightforward. (On the contrary in the 
gauge A + = the amplitude for hard gluon exchange is infrared divergent. The divergent contribution is, however, 
exactly canceled out with that in the fermion propagator in this gauge, cf. Q.) The chosen gauge is convenient to 
evaluate the high momentum component of the pion wave function. For the evaluation of the parton distribution 
within a nucleon another gauge is more appropriate. This mismatch does not introduce, however, additional problems 
because the factorization theorem justifies the possibility to choose independently the gauges for the amplitude of 
pion fragmentation into jets and for the parton distribution within a nucleon. The perturbative tail is obtained as the 
result of the one gluon exchange interaction acting on the soft part of the momentum space wave function, defined as 



ip(k,y) = (h,y I 



111- 



(45) 



By definition, ip is dominated by its non-perturbative low-momentum components. However, the amplitude we 
compute depends on the high momentum tail, x- For this component, perturbation theory is applicable and we use 
the one-gluon exchange approximation to the exact qq wave function of Eq. ( f45| ) to obtain %, valid for large values of 
Kf , as 



X-{k tl x) = -~4irC F 



1 



kf+mj, 



x{l—x) 



dy 



2(2tt) : 



■V°(kt,x;l t ,vMlt,y) 



(46) 



with 



,r gn , \ u{x,k t ) u{y,l t )v{l-x,-k t ) v(l-y,-l t ) 
V y (kt,x;l t ,y) = a s ^— 7m 7= 7== — Ju s= — « M 



0(y - x) 



x l—y 



(kt-lt) 2 
y-x 



- + (x -i 



x,y^l-y) 



(47) 



where Cf 



and d^ v is the projection operator of the gluon propagator evaluated in the light cone cone 



gauge defined above : d^ v = S^u 



?vO)kv+fc M pLO) 

(p;fc) 



Here k is the gluon four-momentum. The range of integration 



over l t is restricted by the non-perturbative pion wave function ip. 

Then in the evaluation of V s we set m q and It to everywhere in the spinors and energy denominators. This is 
legitimate because of a lack of infrared divergences in the amplitude of hard process Jl7j . 

Thus 



V 9 (k u x;l t ,y) 



a s {k 2 t ] 



x(l - .i)y(l - y) 



V(x,y) 



(48) 



2 An early version of this subsection which did not include discussion of renormalization effects has appeared in [ |l7| . 
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where in the lowest order over coupling constant a s but keeping leading power of fcf : 

V(x, y) = 2 [(% - x)x{\ -y) + e(y-x)) + (x^l-x,y^l- y)] , (49) 

We want to draw attention that in the gauge chosen in the paper dominant contribution arises from the components 
of gluon propagator transverse to the directions of pion and nucleon four momenta. 
The net result for the high k t component of the pion wave function is then 

where nonperturbative support 



The analysis of experimental data for virtual Compton scattering and the pion form factor performed in 
shows that this amplitude is not far from the asymptotic one Jl7) for k 2 > 2 — 3 GeV 2 , 

(\>(k\ -> oo,y) = a y(l - y), (52) 

where a = y^/V with /„• w 93 MeV. 

Equation ( |50| ) represents the high relative momentum part of the pion wave function in the lowest order of pQCD 
when running of the coupling constant is neglected. Using the asymptotic function (|5^) in Eq. ( |50| ) leads to an 
expression for \(kt,x) cx a;(l — x)/k 2 

4irC F a s (k?) 

X{k t ,x) = a x(l-x). (53) 

The next step is to use the renormalization invariance of theory to include the k t dependence of the coupling 
constant : 

4-7T 

a s {k 2 t ) = w , (54) 

with /3 = 11 — §«■/■ This can be easily done similar to (4^] where the relationship between (f>(Q 2 ,x) and x(kt,x) and 
QCD evolution equation for (j)(Q 2 , x) have been deduced. 

The quark distribution function 4>{x, Q 2 ) -the amplitude for finding constituents with longitudinal momentum x in 
the pion which are collinear up to the scale Q 2 is 



^ Q ^dM)k x{ktx) ^ (55) 

The factor : , j-^ = (In -tS — ) 1F /& arises from vertex and self-energy corrections. (By definition the running 

a f\ K t> QCD 

coupling constant includes the renormalization factor of the fermion propagator. Such a renormalization factor is 
absent in the definition of x(fc t ,x)). One of advantages of the gauge chosen in this paper is that 7^ 

If = C F , (56) 

has no infrared divergences in difference from |f(| where gauge A + has been chosen. Calculations are significantly 
simplified because the dominant contribution is given by transverse gluon polarizations only. The evolution equation 
for cj)(x,Q 2 ) gives 

00 /in O 2 \ ~~ 1n 

(j>(x,Q)=x(l-x)J2 a nCn(l-2x) ( JJ^~ ) , (57) 

n=0 \ QCD J 

where j n = ^(1 — ( w+1 ) 2 (» +2 ) + ^J2k=2 +1 ^-/^)' — 2x) are Gegenbauer polynomials . Coefficients a n are the 

subject of discussions in the literature and they can be estimated within the current models. The above equation 
makes it possible to calculate the high k t behavior of x(fct, x) as 
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(58) 



Thus the z, n t dependence of the cross section for the diffractive dijet production should be sensitive at moderately 
large k t to the terms involving Gegenbauer polynomials of the order greater than in the pion wave function if a n 
are sufficiently large. The process of the photoproduction of jets at HERA will be appropriate for this purpose. 
It follows from the above equation that the asymptotic wave function is as follows: 

X(h,x) = 47rCj ," s ^^ a x(l - x) (in-jS — J . (59) 
K t \ a qcd ) 

Thus QCD predicts the dependence of x(kti x ) to be of the form used in Ref. ||. 



B. Interaction with the target 

To compute the amplitude Ti it is necessary to specify the scattering operator /. We will fix the transverse recoil 
momentum of the target at zero to simplify the discussion. The transverse distance operator b = (b q — bq) is canonically 
conjugate to Kf At sufficiently small values of b, the leading twist effect and the dominant term at large s arises from 
diagrams in which the pion fragments into two jets as a result of interactions with the two-gluon component of the 
gluon field of a target, see Fig. 1. The perturbative QCD determination of this interaction, which is a type of QCD 
factorization theorem, involves a diagram similar to the gluon fusion contribution to the nucleon sea-quark content 
observed in deep inelastic scattering. One calculates the box diagram for large values of n t using the wave function 
of the pion instead of the vertex for 7* — > qq. The a pp lication of the technology leading to the QCD factorization 
theorem in the impact parameter space leads 18 5(if| to 



n 2 



f(b 2 )=is—b 2 [x N G N (x N ,Q 2 cS ) + 2/3x N S N (x N ,Q 2 cS )] a s (Q 2 cS ), (60) 

in which Gat is the gluon distribution function of the nucleon, Sn is the sea quark distribution function of the nucleus 
for a flavor coinciding with that of the qq dipole, and Q 2 S — A. The factor 2/3 appearing in the second term is 
the same as in the LO approximation for the longitudinal structure function and exclusive vector meson production 
| pl| . The only difference is that in our case the number of flavors is unity. For our kinematics, it is reasonable 
to use X(x = 10~ 3 ) = 9 p5| . The formula (|6(]) should be modified when applied to hard diffractive processes. 
The mass difference between the pion and the final two-jet state requires that the reaction proceeds by a non-zero 
momentum transfer to the target. This means that the function Gn should be replaced by the skewed (or off-diagonal 
or generalized) gluon distribution. Thus the distribution function should depend on the plus components x±,X2 of 
the momenta hi, &2 of the two exchanged gluons, Eq. (R). 

The difference between the skewed and ordinary gluon distribution is calculable in QCD using the evolution equation 
for the skewed parton distributions pT]j52| . The kinematical relation between x\ and X\ — x-i is given in Eq. (^|). But 
x\ is close to xn of Eq. (|6l|), while X2 is small in the calculation of T\. The skewed parton distribution can be 
approximated by a gluon distribution |53|]54|] if 



X ^^ +X2) I^M&Ws- (61) 

While including the effect of skewedness would alter any detailed numerical results, the qualitative features of the 
present analysis would not be changed. 

The most important effect shown in Eq. (|6^) is the b 2 dependence, which shows the diminishing strength of the 
interaction for small values of b. To simplify formulas it is convenient to rewrite a in the form: 

Kb 2 ) = ™^b 2 = lS ^L (-V 2 K ) , (62) 

in which the logarithmic dependence of a s and the gluon distribution on b 2 are included in 0$. It is easy to check 
by direct calculations of Feynman diagrams that the operator V 2 acts on the transverse momentum variables of the 
pion wave function. Our notation is that (6q) represents the pionic average of the square of the transverse separation, 
and within the leading log accuracy 
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°0 ~ C 2\r ^(skewed)/ 2\i 

» ya,(K t )[xjv Gjv (^l,^,^)], (63) 

in which the ordinary gluon distribution is used as the initial condition for the QCD evolution equation for the 
skewed/generalized gluon density. This result is a factor of four smaller than presented in Ref. fa]. 



The result (63) holds for xm about 10 -2 . For xm of about 10~ 3 or smaller, the second kinematic regime mentioned 
in the introduction is relevant, and one would obtain different results. For still smaller values of x, say x ~ 1CP 5 
non-linear gluonic effects become important, and the present treatment of the qq interaction with the target may be 
insufficient. 

C. One Gluon Exchange in the Pion— Ti 

The necessary inputs to evaluating T\ are now available. The approximate pion wave function, valid for large 



relative momenta, is given by Eq. (J53J). The interaction / is given by Eq. (62). The use of Eq. (62) allows a simple 
evaluation of the scattering amplitude T\ because the b 2 operator acts on the pion wave function (here ero is treated 
as a constant) as — V 2 t , leading to the result 

a 4irC F a s {n 2 t ) n k 2 t Qf 
Ti = -Ais-—- j ( ln T2 )' a z(l-z). (64) 



(b 2 ) n\ A? 



QCD 



This amplitude T\ is of the same form as our 1993 result ||. The present result is obtained directly from QCD, in 
contrast with the earlier work which used some phenomenology for the pion wave function. 

Corrections to Eq.(|6j) are of the order — For example, a literal application of Eq. ( |62| ) would lead to a factor 



In ■ 



'-t . However, similar corrections may arise from other effects not considered here. So a calculation of such 



corrections is beyond the scope of this paper. 

/ k 2 \2C F /6eta 

Our Kt dependence of T\ leads to: da ^ K ^ cx — QCD ./ for xn ~ 10 -2 . This can be understood using 

simple reasoning. The probability of finding a pion at b < is oc b 2 , while the square of the total cross section for 
small-dipole-nucleon interactions is oc b . Hence the cross section of productions of jets with sufficiently large values 

2C F 

Qs (fc?) 2 (l„(^)^ 

of n t integrated over d 2 k t is oc — K a QCD , for xn ~ 10 2 , leading to a differential cross section oc This 

reasoning ignores the dependence of the gluon structure function on Kt. For sufficiently small values of x(x ~ 10~ 3 ), 
gluon evolution would give a somewhat different behavior. 

V. NUCLEAR DEPENDENCE OF THE AMPLITUDE 

The picture we have obtained is that the amplitude is dominated by a process in which the pion becomes a qq pair 
of essentially zero transverse extent well before hitting the nuclear target. This point-like configuration (PLC) can 
move through the entire nucleus without expanding. The qq can interact with one nucleon and can pass undisturbed 
through any other nucleon. For zero momentum transfer q t to the nucleus, the amplitude M. (A) takes the form 



M{A) = A M{N) A ^ K _ f J J J> 1 + 2 AV3 ) = AA^(iV)r , (65) 



G A (xi,x 2 ,m 2 f ) 
' AGN(xi,X2,m 2 ) V" <b 2 >K 2 

in which the skewedness of the gluonic distribution is made explicit, and where the real number e > 0. Observe the 
factor A which is the dominant effect here. This factor is contained in the ratio of gluon distributions in a nucleus 
and in a nucleon Q. This dependence on the atomic number is a reliable prediction of QCD in the limit m 2 and 
s — > oo, with fixed xn (of Eq. (|6l"|)). (The quantity rrif is the mass of the dijet system.) On the contrary, for x^ — > 
with fixed to 2 , the nuclear shadowing of the gluon distribution becomes very important ]5],|36|. 

The e correction term in Eq. ( p5| ) is a higher twist contribution which arises from a single rescattering which can 
occur as the PLC moves through the nuclear length (Ra oc A 1 / 3 ). That e > was a somewhat surprising feature of our 
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1993 calculation because the usual second order rescattering, as treated in the Glauber theory, always reduces cross 
sections. This highly unusual sign follows from the feature in QCD that the relative contribution of the rescattering 
term (screening term) decreases with increasing size of the spatially small dipole. The key features of the usual first 
order term are / = ia , and those of the usual second order term are if 2 = —ia 2 . Note the opposite signs. For us 
here / = ia$b 2 / < b 2 >, For very large values of k 2 the operator b 2 , as applied to the pion wave function Xi which 

falls with k 2 as gives b 2 x = — So the first term fx — — «4cro/(< b 2 > k 2 )x now has the same sign as the 

K t 

second-order term: if 2 x = —i[<^o/ < b 2 > b 2 ] 2 x = ~ i32[<7o/(< b 2 > n t )] 2 X 
The differential cross section takes the form 

dq 2 dt ' 1 j 

for small values of t. Note that 

-t = q 2 - < min , (67) 
where — i m i n is the minimum value of the square of the longitudinal momentum transfer: 

/ m 2 f - m 2 \ 2 

Our discussion below is applicable in the kinematics where —t m i n R 2 A /?> < 1 so that the entire dependence of the cross 
section on t m i n is contained in the factor e'min-fiU/ 3 . 
One measures the integral 

a(A) = J dt^ = ^rA 2 T 2 a(N) . (69) 
A typical procedure is to parametrize a (A) as 

a(A) = G\A a (70) 

in which a\ is a constant independent of A. For the Ra corresponding to the two targets Pt (A = 195) and C(A = 
12) of E791, one finds a ~ 1.45. The experiment |l4j does not directly measure the coherent nuclear scattering cross 
section. This must be extracted from a measurement which includes the effects of nuclear excitation. The extraction 
is discussed below. 

As pointed out previously p7[ , the values of our multiple scattering correction of our 1993 calculation ||] were 
overestimated by a factor of approximately four. This is because the Co was chosen to be larger by a factor of 4 
than in |Q and in Eq. (J60|) . We now find that for values of n t greater than about 2 GeV, the coefficient a could be 
enhanced by between 0.0 and 0.08, depending on the value of Kt- Taking 0.04 as a mean value one finds a w 1.5. 
This estimate depends on the use of a model for the non-perturbative part of \ir), and also on the validity of a simple 
eikonal treatment for the multiple-scattering corrections which is questionable at the high energies that we consider 
here. 

Another potentially important A-dependent effect is the nuclear shadowing of the parton densities. Very little direct 
experimental information is available for the A-dependence of the gluon structure function. The analyses of the data 
combined with the momentum sum rules and the calculation of gluon shadowing at small x suggest that the value 



x ~ 0.01 (which corresponds to the kinematics of 14|) is in a transition region between the regime of an enhancement 



of the gluon distribution at x ~ .1 and the strong shadowing at x < 0.005 |20|-|23j|. Using the A-dependence of F2A 
as a guide, and in particular the NMC ratio F2s n /F 2 c the shadowing may reduce a for the Q kinematics by 
Aa 0.08. 



VI. EXPERIMENTAL CONSIDERATIONS 

The requirements for observing the influence of color transparency were discussed in 1993 The two jets should 
have total transverse momentum to be very small. The relative transverse momentum should be greater than about 
> 2 GeV and the mass of the diffractive state should be described by the formula 
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Maintaining this condition is necessary to suppress diffraction into a qqg pair in which the gluon transverse momentum 
is not too small. 

It would be nice if one could measure the jet momenta precisely enough so as to be able to identify the final 
nucleus as the target ground state. While it is very feasible to consider for eA colliders, this is impossible for high 
energy fixed target experiments, so another technique must be used here. The technique used in is to isolate the 
dependence of the elastic diffractive peak on the momentum transfer to the target, t, as the distinctive property of 
the coherent processes. This was done by first introducing a cut on the momentum of the observed pions by requiring 
that they carry more than 90% of the incident momentum. This sample was than analyzed as a function of the total 
transverse momentum of the system. A strong coherent peak was observed with the slope consistent with the coherent 
contribution. The background was fitted as a sum of the coherent peak, inelastic diffraction with a nuclear break up 
and the term due to inelastic events where some hadrons were not detected. 

The amplitude for the non-spin flip excitations of low-lying even-parity nuclear levels ~ — t, due to the orthogonality 
of the ground and excited state nuclear wave functions. Thus the cross section of these kinds of soft nuclear excitations 
integrated over t is suppressed by an additional factor of l/R\ ~ A~ 4 / 3 compared to the nuclear coherent process. For 
y/—i Ra 3> 1 where q is the momentum transfer to the nucleus q = \f~-i the background processes involving nuclear 
excitations vary as A, so an unwanted counting of such would actually weaken the signal we seek. For qRA 3> 1 the 
diffractive peak cannot be used as signature of diffractive processes to distinguish them from non-diffractive processes 
whose cross section oc a(irN) oc ^4 75 . Thus substantial A-dependence, a cx A a (with a « 1.5), as predicted by QCD 
for large enough values of n tl should be distinguishable from the background processes. 

The amplitude varies as M.(A) ~ a s / 1 k\ 

M{A) ~ a s x N G N (x N , Ql ff )/4, (72) 

where Q 2 ff ~ ^ K t- For the kinematics of the E791 experiment, where increases oc ftf, the factor 
a s XffG n(%n iQef f) ^ s a rather weak function of K t . For example, if we use the standard CTEQ5M parameteri- 
zation we find <t(A) ~ l/«f' 5 for 1.5 < n t < 2.5 GeV which is consistent with the data |L4]]. For the amplitude 
discussed here, a (A) ~ (z(l — z)) 2 which is in the excellent agreement with the data 0. 

A. Extracting the coherent contribution 

The experiment is discussed in Ref. fl4{| . The main advantage of this experiment is the excellent resolution of 
the transverse momentum. The reference also shows the identification of the dijet using the Jade algorithm, and it 
displays the identification of the diffractive peak by the q 2 dependence for very low q 2 . This dependence is consistent 
with that obtained from the previously measured radii R{C) = 2.44 fm, and Rp t — 5.27 feu. The key feature is the 
identification of the coherent contribution from its rapid falloff with t. 
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FIG. 13. Contributions to the total nuclear diffractive cross section. The wavy lines in this figure denote amplitude for the 
scattering of qq pair by a nucleon. (a) The terms with i = j. (b) i 7^ j. 



We discuss the extraction of this signal in some detail. We consider the contribution to the total cross section 
that arises from the diffractive production of the dijet. The total cross section includes terms in which the final 
nucleus is not the ground state. The nuclear excitation energy is small compared to the energy of the beam, so that 
one may use closure to treat the sum over nuclear excited states. Then the cross section is evaluated as a ground state 
matrix element of an operator Yli,j e ^ '( Ti ~ Tj ' > = A + Ylijtj e l( l'( r< ~ r " ; see Fig. 13. The result, obtained by using Xi 
relative to the nuclear center of mass, and by neglecting correlations in the nuclear wave function is given by 



do a 
dt 



A + A(A - l)F 2 A ( t 1 



A-l 



The factor in the argument of Fa is due to accounting for nuclear recoil in the mean field approximation, cf. 
|p5| . This formula should be very accurate, for the small values of t relevant here. The contribution of the coherent 
processes to the total cross section is given by 

J —coherent -V_ 

and the contribution of excited nuclear states is the difference: ^f- — da dt , which vanishes at t = 0. The 

experiment proceeds by removing a term oc A from ^f- which has no rapid variation with t. This defines a new cross 
section which is actually measured experimentally. 

t-^-M*^)^ (75) 

The integral of this term over t can be extracted from the data: 

- f d °A __ ?>A{A - 1) da N 3(A - l) 2 da N 

ai ~J dt rjf+R^ dt | t= o W r%+R\ dt | i=0 ' [<b) 



Here the factor rj^ takes into account the slope of the elementary cross section, assuming that it is determined solely by 



the nucleon vertex. Note that the result (|76|) differs by a factor of ^ A A P from the A-dependence predicted previously 
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for coherent processes, recall Eq. (|66|). Using A— 12, 195, the nuclear radii mentioned above, rjv = 0.83/m and fitting 
the ratio of cross sections obtained from Eq. (|76| ) with the parametrization a oc A a , gives then 



instead of a — 1.45. 

The result |14j of the experiment is 



a = 1.54, (77) 



a w 1.55 ± 0.05 , (78) 



which is remarkably close to the theoretical value shown in Eq. (77). The sizes of our multiple scattering and nuclear 
shadowing corrections, which work in the opposite directions, and which were discussed in the previous Section, are 
of the order of the experimental error bar. 

VII. ELECTROMAGNETIC BACKGROUND 

Because very low values of q 2 are involved, one could ask if the process occurs by a one-photon exchange (a type 
of Primakoff effect) instead of a two-gluon exchange. If the momentum transfer is very low, the process is peripheral 
and there would be no initial or final state interactions. Thus, it is necessary to estimate the relative importance of 
the two effects. 

This amplitude is caused by the exchange of a virtual photon of four-momentum q (q 2 = t) with the target. The 
nuclear Primakoff amplitude is then given by 

M P {A) = e^^M (PA + Pir^F A (t) « 2eMJ- ■ ^|«>^, (79) 

where we will use the decomposition —t = q 2 =q 2 — t m - m . A photon can be attached to any charged particle, so a 
direct calculation would involve a complicated sum of diagrams. We may simplify the calculation by using the fact 
that the electromagnetic current is conserved. This application is simplified by the use of Sudakov variables, which is 
a necessary first step. Accordingly, we write 

q = + j3p w + q t . (80) 

Conservation of four-momentum gives 



q 2 -mj 



^W^pZY " . ■ (81) 



Then conservation of current can be written as 



<tt| J om • q\qq) a a<7r| J cm • ^\qq) + J cm • p«\qq) - (tt| J cm • q t \qq) = 0. (82) 



The use of Eq. ( plf ) and keeping only the leading term in /i 2 /s, where \x is the typical mass involved in the considered 
process, allows us to neglect the (3 term of Eq. ( |82"| ) so that 

a(7r\J™ ■ ^\qq) = (ir\J™ ■ q t \qq) (83) 

By definition, the transverse momentum of a pion is zero, so the dominant (in terms of powers of Kt) contribution in 
Eq. (p3] ) is given by photon attachments to quark lines, and the matrix element is given by ^ 

<vr| J cm • q t \qq) = K t )q t ■ K t (2/3z - 1/3(1 - z)). (84) 



3 The above ( |84[ ) differs from that of the first version of our paper which appeared in hep-ph. We are indebted to D.Ivanov 
and L.Sczymanowsky - who drew our attention to the misprint in this version of the paper. 
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The generalization of this result to account for all Feynman diagrams having the same powers of s and k 2 is almost 
trivial. The relative contribution of other diagrams is oc k' t /nt where k' t is the transverse momentum of quarks in the 
intermediate state. But within the a s In k 2 /Aqcd approximation k' 2 <C k 2 so this contribution does not lead to a 
\u.k 2 /Aq CD term. Thus the above formula is valid within the a s In k 2 /Aq CD approximation when q s « 1. 
The net result, obtained by using Eq. (|83|) in Eq. (f79|), is 



M P {A) = e2 f^' Kt)Z F A (t)^2q t ■ n t (2/3 - z) (85) 

Qt ^min f^t 

which should be compared with the amplitude of Eq. (Q) (including the effect of the nuclear form factor, which enters 
at non-zero values of t, but ignoring the logarithmic correction) written as 

M(A) = -i X7T (z, Kt )A^^-F A (t) « 4ixv(z, n t )/k 2 t A-^-F A {t) (86) 
< b z > < b z > 

The Primakov term has been evaluated also in the paper of D.Ivanov and L.Sczymanowsky p2| ] . As the consequence 
of different approximations result obtained in J62[ differs from ours. We use conservation of the e.m. current to 
deduce the Weizecker- Williams approximation and to account for the cancellation between diagrams corresponding 
to attachments of photon to the different charged constituents of a pion. After accounting for the cancellation we 
restrict ourselves by the contribution of diagrams enhanced by the large factor n t from the vertex of quark(antiquark)- 
photon interaction. Dominance of photon interaction with external quark lines in the pion fragmentation into 2 jets is 
another form of the QCD factorization theorem which properly accounts for the conservation of the e.m. current and 
the gauge invariance of QCD. On the contrary (6^] interacting particles are put on mass shell before the separation of 
scales and the cancellations between different photon attachments has been taken into account. These approximations 
have problems with the conservation of the e.m. current and the renormalization group in QCD. 

The ratio of electromagnetic and strong amplitudes is given by 

M P {A) _ , 2 Z 2/3 -z g.K t 

M(A) A a / <b 2 > 2(q 2 - i min ) {0 '> 

Using gf fa 0.02 GcV 2 (the smallest value measured in Z/A=l/2, e 2 = 47r/137,K t = 2 GeV, a / < b 2 >ps 2.5 
(this is 1/4 of the value of Rcf. ||), and taking q t parallel to n t we find that 

« -1/8(2/3 - z) i « -0.02 i, (88) 

with z = 1/2. 

Thus the Primakoff term is very small and, because of its real nature, interference with the larger strong amplitude 
(which is almost purely imaginary) is additionally suppressed. We may safely ignore this effect for the energy range of 
Rcf. or any contemplated fixed target experiment. At collider energies, it will be possible to measure jets of much 
larger values of transverse momentum, so any complete theoretical analysis should account for this electromagnetic 
interaction. 

For a heavy nuclei target another electromagnetic process oc Z 2 a 2 m in the amplitude gives a contribution. This is 
a dijet production due to the two-photon exchange, a version of Fig. 8 in which the exchanged gluons are replaced 
by exchanged photons. For small transverse momenta of quarks I 2 -C k 2 in the pion wave function, this contribution 
is suppressed as the power of s in the amplitude. This is because, in the calculation of the imaginary part of the 
diagram, x-i for the exchanged photon is given by 

k 2 

x 2 v=--^> (89) 

except for a very narrow region of z near z = 0, 1, which is suppressed by the decrease of the pion wave function. Thus 
X2 < and according to our previous arguments, this contribution should be very small. For large It this contribution 
may also be expressed in terms of the same pion wave function as in the case of the two-gluon exchange. It is easy to 
estimate this EM amplitude obtained from the two-photon exchange: 

M(rr) = a 2 em {s/k 2 ) x {z, k t )„Z 2 J f^F A (l t )F A ( qt - l t )(8e iea )(l - 2z). (90) 

Here ei(e2) is electric charge of quark (anti-quark) in the units of the electric charge of electron. The lower limit of 
integration over I 2 is (k 2 /2z(1 — z)) 2 Thus the contribution of this term to the cross section should have the same 
z and k 2 dependence as the two-gluon exchange term, but with a much faster Z dependence (oc Z 4 ). This term is 
negligible also. 
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VIII. DISCUSSION AND SUMMARY 



The use of the experimentally measured |Q value of a — 1.55 (recall Eq. (|70|)) leads to 

a(Pt) 



a(C) 



75 . (91) 



The typical usual nuclear dependence of the soft diffractive processes observed in the high energy processes is « A 2 / 3 
or a = 2/3. The use of a Glauber approximation with a typical hadronic cross section for the final system tends to 
predict the A dependence as « A 1 / 3 . An account of color fluctuations J56| predicts « A 2 / 3 , in agreement with the 
FNAL data |57]] which would give 

ousualCP^ = 7 (92) 

CUSUAl(O) 

Thus color transparency causes a factor of 10 enhancement! This seems to be the huge effect of color transparency 
that many of us have been ho ping to find. It is also true that, as noted in the Introduction, that the n t and z 
dependence of the cross section |14| is in accord with our prediction. 

All of this looks very good, but it is necessary to provide some words of caution. Our analysis was related to 
a nuclear coherent process involving a qq final state. If the experimental signal is significantly contaminated by 
incoherent nuclear effects or by qqg final states, our analysis might not be applicable. However, the experimental |14| ] 
extraction of the coherent peak using the q 2 dependence of the amplitude, and the measurement of the two-jet (as 
opposed to three-jet) cross section seem very secure to us, except for the small correction discussed in Sect. V. Another 
worry is that the color transparency effect seen in Ref. jl4| seems to start for values of n t near 1 GeV. These are lower 
than suggested in Ref. ||. These earlier predictions used modeling of non-perturbative effects, and such modeling 
may be necessary to guess the lowest values of Kt for which color transparency would occur. The reasoning of the 
present paper uses perturbative QCD, which becomes more reliable as n t increases. This is because the competing 

amplitudes T 2 3 4 are decreased relative to Ti by a factor of ^ (~ -04 for n t — 1 GeV) or a s (n 2 ). A coherent sum 

1 ' K t 

of the sub-dominant amplitudes could provide a significant correction to our dominant pure amplitude. However, the 
observed falloff of the cross section with K t , combined with the z and A dependence, does provide very strong evidence 
for color transparency. 

It is worth noting similarities and differences between the process we discuss in this paper and another factorizable 
process, that of hard exclusive electroproduction of mesons |3£],|5^]. Both processes allow a simple geometric inter- 
pretation in the transverse coordinate representation: a convolution of the initial wave short-distance wave function, 
^in(-2,6), the dipole-target cross section, a(b,x), and the final wave function ipfin(b, z). However, in the case of the 
vector meson production, the Q 2 dependence of the longitudinal photon wave function at small b causes the final 
vector meson wave function to be evaluated at b 2/Q. In the case of the pion dijet diffraction the pion wave 
function enters at small values of b ~ — . 

If color transparency has been correctly observed in the n + A — > qq + A (ground state), there would be many 
implications. The spectacular enhancement of the cross section would be a novel effect. The point-like configurations 
would be proved to exist. This would be one more verification of the concept and implications of the idea of color. 
Furthermore, the definitive proof of the existence of color transparency means that we now have available a new 
effective tool to investigate microscopic hadron, nuclear structure at hundreds of GeV energy range. At lower energies 
of a ten's of GeV color transparency is masked to some extent by the diffusion of a spatiallysmall quark-gluon wave 
package to the normal hadronic size vM- Still it is possible that previous experiments |59[] showing hints pOfl of 



color transparency (for a review see Ref 12 1) probably do show color transparency. Efforts j6TJ to observe color 
transparency at Jefferson Laboratory, and at HERMES (DESY) should be increased. The electron-ion collider would 
provide numerous possibilities for studying color transparency both in di-, tri-jet coherent production as well as in 
exclusive processes. 

The observation of CT confirms the idea that the life span of the perturbative phase can be increased by the large 
Lorentz factor associated with high energy beams. A challenging problem would be to explore this idea to observe the 
perturbative phase in a "macroscopic" volume. One manifestation of this would be the production of huge blob-like 
configurations of Huskyons |i3| . These different configurations have wildly different interactions with a nucleus |34| , 
so that the nucleon in the nucleus can be very different from a free nucleon. More generally, the idea that a nucleon 
is a composite object is emphasized by these findings. Some configurations of the nucleon interact very strongly with 
the surrounding nucleons; some interact very weakly. This means that the nucleon in the nucleus can be very different 
from a free nucleon. This leads to an entirely new view of the nucleus, one in which the nucleus is made out of 
oscillating, pulsating, vibrating, color singlet, composite objects. 
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The technical purpose of this paper has been to show how to apply leading-order perturbative QCD to computing 
the scattering amplitude for the coherent processes: nN — > J J N and nA — > J J A. The high momentum component 
of the pion wave function, computable in perturbation theory, is an essential element of the amplitude. The dominance 
of the amplitude of the T\ term of Eq. (64) is obtained by showing that the corrections to it, which at first glance seem 
to be of the same order in the coupling constant, are vanishingly small. This vanishing, obtained using arguments 
based on analyticity, causality, and current conservation, is equivalent to the verification of a specific space-time 
description of the event: the pion produces its point-like component at distances well before the target. Furthermore, 
for the conditions of the experiment fl4| studied here, the competing electromagnetic production process is shown 
to yield a negligible effect. It seems that perturbative QCD can be applied to the coherent nuclear production of 
high-relative momentum dijets by high energy pions. 

It therefore seems interesting to consider similar reactions involving other projectiles such as photons, kaons, and 
protons. The observations of the coherent photoproduction of the J/?/; from nuclear targets has long been known |65| ] 
to have an A dependence which is very similar to that observed here, but the authors of [ |65| did not interpret it as 
color transparency phenomenon. Later on HI and ZEUS detectors at HERA investigated exclusive photoproduction 
of the J/ip meson from a proton target. The striking qualitative predictions for this process based on the QCD 
factorization theorem, such as energy and t dependence, are in accordance with the data; for a recent review see 
[0,^6|. Thus now there exist serious reasons to believe that color transparency phenomenon has been observed in 
the combination of coherent photoproduction of J/ip from nuclei(FNAL) and from the nucleon (HERA). Our present 
theory can be used for kaon projectiles with little modification. Because the kaon has a smaller size than the pion, 
we expect that the amplitude for a kaon-induced process should be somewhat larger than that of the pion induced 
process discussed here. The same analysis should be applicable to the photoproduction of high n t qq pair of light 
quarks. It is important here that the contribution of bare photon coupling to the qq is proportional to quark's bare 
mass and is negligible for forward scattering |31p^ ] . The contribution of a target gluon with kf t ss k\ discussed in |3|] 
is suppressed by Sudakov and w-i form factors discussed above. However, the PQCD physics of light quarks will be 
masked to some extent by another striking prediction of QCD, which is the enhancement of the diffractive production 
of charmed dijets because of a large bare mass of charmed quarks. 

The study of high energy coherent production of jet systems from nuclear targets seems to be a very productive 
way to investigate both perturbative QCD and microscopic nuclear structure by exploring the diverse effects of color 
transparency. Such studies in the region of hundreds of GeV seem ideally suited for the non-destructive investigation of 
a microscopic hadron, and nuclear structure. It may be possible to remove a piece of hadron (qq pair..) or to implant 
(strangeness in the center of a nucleus...) without destruction of a target. Such investigations resemble modern 
methods of surgery which avoid cutting muscles. So it seems appropriate to name this new field of investigations as 
micro-surgery of a hadron, or of a nucleus. 
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